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ABSTRACT 


J 

A  mathematical  model  is  developed  to  determine  the  heat 
transfer  rate  of  a  cylindrical  condenser  section  of  a  rotating 
heat  pipe.  This  model  is  coupled  to  an  existing  code  and  an 
analysis  is  accomplished  on  both  a  smooth  and  axially  finned 
condenser.  The  results  of  this  analysis  are  compared  to  those 
of  a  similar  analysis  performed  on  a  tapered  condenser  heat 
pipe  using  identical  geometric  and  operating  parameters. 

Results  of  the  comparison  indicate  cylindrical  condensers 
have  a  lower  heat  transfer  rate  than  equivalent  tapered  con¬ 
densers.  This  reduction  in  heat  transfer  rate  is  due  to  a 
greater  condensate  film  thickness  and  is  most  significant  in 
a  smooth  condenser. 

Axially  finned  condensers  with  triangular  and  rectangular 
fin  profiles  are  also  compared.  The  rectangular  fins  are 
assumed  to  have  adiabatic  tips.  Results  indicate  the  heat 
transfer  rates  for  these  two  profiles  vary  by  only  0.40  per 
cent  for  both  tapered  and  cylindrical  condenser. 
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I.  INTRODUCTION 


A.  THE  ROTATING  HEAT  PIPE 

The  rotating,  wickless  heat  pipe  is  a  closed  container 
designed  to  transfer  large  amounts  of  heat  from  rotating 
machinery  components.  Essentially,  it  consists  of  three 
main  components:  a  cylindrical  evaporator  section,  a  con¬ 
denser  section  which  may  be  either  tapered  or  cylindrical  in 
shape,  and  a  fixed  amount  of  working  fluid.  A  typical  tapered 
rotating  heat  pipe  is  shown  in  Figure  1. 

When  the  heat  pipe  is  rotated  about  its  longitudinal  axis 
at  a  speed  above  a  certain  critical  value,  the  working  fluid 
forms  an  annulus  in  the  evaporator  section.  Note  in  Figure  1 
that  the  diameter  of  the  evaporator  is  larger  than  the  con¬ 
denser.  This  larger  diameter  provides  a  greater  liquid  reser¬ 
voir.  As  heat  is  added  to  the  evaporator,  the  fluid  in  the 
evaporator  will  vaporize.  The  vapor  will  flow  axially  towards 
the  condenser  as  a  result  of  a  slight  pressure  difference, 
transporting  the  latent  heat  of  vaporization  with  it.  In  the 
condenser  end,  external  cooling  of  the  condenser  causes  the 
vapor  to  condense.  In  the  case  of  a  tapered  heat  pipe,  the 
centrifugal  force  due  to  the  rotation  of  the  pipe  has  a  com¬ 
ponent  acting  along  the  condenser  wall  which  accelerates  the 
liquid  condensate  back  to  the  evaporator  to  complete  the  cycle. 


condensate 


A  cylindrical  heat  pipe,  on  the  other  hand  relies  on  a  hydro¬ 
static  pressure  gradient  to  drive  the  liquid  condensate  back 
to  the  evaporator. 

B .  BACKGROUND 

The  first  theoretical  investigation  into  the  performance 
of  a  tapered  rotating  heat  pipe  at  the  Naval  Postgraduate 
School  was  accomplished  by  Ballback  [Ref.  1]  in  1969.  He 
examined  the  limits  in  heat  transfer  controlled  primarily  by 
fluid  dynamic  considerations.  In  particular,  he  considered 
the  following  four  limits  on  heat  pipe  performance:  a)  the 
boiling  limit,  b)  the  entrainment  limit,  c)  the  sonic  limit 
and  d)  the  condensing  limit.  Tantrakul  [Ref.  2]  calculated 
these  limits  for  a  specific  heat  pipe.  He  found  the  condens¬ 
ing  limit  was  the  controlling  limitation.  In  fact,  the  cal¬ 
culated  heat  transfer  rate,  based  on  the  condensing  limit  was 
l/10th  the  heat  transfer  rate  for  the  next  lowest  limit,  the 
entrainment  limit. 

In  order  to  overcome  this  condensing  limitation  and  thus 
increase  the  heat  transfer  rate  of  the  rotating  heat  pipe, 
the  concept  of  an  internally  finned  tapered  rotating  heat 
pipe  was  considered  by  Schafer  [Ref.  3],  Schafer  developed 
an  analytical  model  for  this  tapered  heat  pipe  with  a  trian¬ 
gular  fin  profile  as  shown  in  Figure  2.  He  assumed  one  dimen¬ 
sional  heat  conduction  through  the  wall  and  fin.  Corley  [Ref. 
4]  developed  a  two-dimensional  heat  conduction  model  using  a 
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Finite  Element  Method  formulation  for  the  same  geometry.  To 
overcome  the  problem  of  few  nodal  points  along  the  fin  surface, 
Corley  assumed  a  parabolic  temperature  profile  along  the  fin 
surface.  Tantrakul  [Ref.  2]  modified  Corley's  computer  code  by 
increasing  the  number  of  finite  elements  from  two  to  three  in 
order  to  minimize  the  error  at  the  apex  of  the  fin.  Purnomo 
[Ref.  5]  developed  a  two-dimensional  Finite  Element  Method 
solution  to  the  steady  state  heat  conduction  problem  using  a 
linear  triangular  finite  element.  Davis  [Ref.  6]  modified 
Purnomo 's  code  to  make  it  compatible  with  COPES/CONMIN  [Ref.  8], 
an  optimization  program.  Davis,  in  his  modification,  found  a 
coding  error  in  Purnomo' s  [Ref.  5]  code,  that  once  corrected, 
permitted  Purnomo 's  corrected  code  to  converge  to  Schafer's 
[Ref.  3]  results. 

Purnomo's  [Ref.  5]  code  is  limited  in  that  it  is  restricted 
to  one  particular  condenser  geometric  configuration,  namely: 
an  axially  finned  tapered  condenser  heat  pipe  with  a  triangular 
fin  profile.  In  that  tapered  finned  condenser  are  difficult 
to  manufacture,  it  is  doubtful  that  any  widespread  practical 
application  of  this  geometric  configuration  will  result.  Cylin¬ 
drical  condensers,  on  the  other  hand,  can  be  manufactured  with 
much  less  difficulty  and  might  find  practical  application. 

This  being  the  case,  a  more  practical  and  beneficial  code 
would  be  one  that  could  analyze  cylindrical  condenser  rotating 
heat  pipes,  both  finned  and  smooth.  In  actuality,  the  most 
beneficial  code  would  be  one  that  could  analyze  the  following 
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four  geometric  configurations:  1)  tapered-internally  finned, 

2)  tapered-smooth,  3)  cylindrical -internally  finned,  and  4) 
cylindrical-smooth.  The  theoretical  heat  transfer  performance 
of  the  four  geometries  could  then  be  compared  to  determine  the 
advantages  and  disadvantages  of  each  design.  An  additional 
advantage  would  be  gained  if  different  fin  profiles,  i.e., 
triangular  vice  rectangular,  could  also  be  analyzed  and 
compared . 

C.  THESIS  OBJECTIVES 

The  objectives  of  this  thesis  are: 

1)  Develop  analytical  models  for  both  cylindrical-smooth 
and  cylindrical-axially  finned  condensers. 

2)  Develop  solution  techniques  to  these  analytical  models 
that  will  account  for  temperature  variations  along  the 
axial  length  of  the  condenser. 

3)  Modify  Purnomo's  [Ref.  5]  code  to  provide  a  solution  to 
the  two-dimensional  steady  state  conduction  heat  transfer 
problem  for  the  following  four  geometric  configurations: 
a)  tapered-smooth,  b)  tapered-finned,  c)  cylindrical- 
smooth,  and  d)  cylindrical - finned . 

4)  Modify  Purnomo's  [Ref.  5]  code  to  provide  the  additional 
capability  of  analyzing  a  rectangular  fin  profile  with 
an  adiabatic  tip. 

5)  Obtain  and  compare  results  of  the  four  geometric  con¬ 
figurations  given  above  for  various  operating  conditions. 
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THEORETICAL  ANALYSIS  FOR  A  CYLINDRICAL  HEAT  PIPE 


A.  INTRODUCTION 

In  a  cylindrical  condenser  heat  pipe,  the  radius  of  the 
condenser  is  constant  along  the  axial  length  of  the  condenser. 
The  flow  of  the  condensate  in  the  absence  of  vapor- liquid 
interfacial  shear,  is  dependent  upon  the  variation  in  hydro¬ 
static  pressure  with  changes  in  film  thickness  along  the  sur¬ 
face  of  the  heat  pipe.  Leppert  and  Nimmo  [Refs.  8  and  9] 
investigated  the  phenomenon  of  film  condensation  on  a  flat 
horizontal  plate.  This  situation  is  similar  to  film  conden¬ 
sation  on  the  inside  surface  of  a  rotating  cylindrical  con¬ 
denser.  In  the  case  of  a  cylindrical  condenser,  the  body 
force,  rather  than  being  the  force  of  gravity,  is  now  the 
centrifugal  force  caused  by  the  rotation  of  the  heat  pipe. 
Weigenseil  [Ref.  10]  and  Tantrakul  [Ref.  2]  compared  experi¬ 
mental  results  for  a  cylindrical  condenser  rotating  heat  pipe 
with  the  theoretical  results  of  Leppert  and  Nimmo  [Refs.  8 
and  9]  and  found  good  agreement.  The  Leppert  and  Nimmo  solu¬ 
tion  was  limited  in  that  it  was  based  on  a  constant  surface 
temperature  along  the  length  of  the  plate.  A  rotating  heat 
pipe,  in  actuality,  has  a  temperature  variation  along  the 
axial  length  of  the  condenser  which  in  some  cases,  may  be 
significant.  This  being  the  case,  it  was  necessary  to  develop 
a  mathematical  model  which  would  consider  the  axial  temperature 


variation  in  the  solution  of  the  heat  transfer  analysis.  In 
the  mathematical  development  that  follows,  a  cylindrical 
smooth  (unfinned)  condenser  will  first  be  considered  in  that 
it  is  the  simplest  case.  The  model  will  then  be  extended  to 
include  a  cylindrical  axially  finned  condenser. 

B.  THEORY  FOR  A  CYLINDRICAL  SMOOTH  CONDENSER 
1 .  Assumptions 

In  developing  the  theoretical  analysis,  the  following 
assumptions  are  made: 

a)  Film  condensation,  not  dropwise  condensation  occurs  in 
the  condenser. 

b)  The  condensate  film  undergoes  laminar  flow. 

c)  Momentum  changes  through  the  condensate  are  small. 

Thus,  there  is  essentially  a  static  balance  of  forces. 

d)  The  vapor  exerts  no  drag  in  the  condensate;  there  is 
no  interfacial  shear. 

e)  The  temperature  distribution  within  the  film  is  linear. 

f)  The  vapor  space  is  essentially  at  one  pressure,  Pv. 

g)  The  density  of  the  fluid  is  much  greater  than  the  den¬ 
sity  of  the  vapor.  Thus,  the  density  of  the  vapor  can 
be  neglected. 

h)  The  centrifugal  force  is  much  greater  than  the  force 
of  gravity  and,  thus,  gravity  may  be  neglected. 


i)  Velocity  gradients  in  the  circumferential  direction 
relative  to  the  pipe  wall  are  negligible. 


j)  The  condensate  film  thickness  is  much  less  than  the 
radius  of  curvature  of  the  condenser  wall. 

k)  The  rotating  heat  pipe  is  operating  at  steady  state 
conditions . 

2 .  Condensate  Momentum  Equation  (X-Direction) 

By  applying  the  above  assumptions  and  the  coordinate 
system  shown  in  Figure  3,  an  analysis  similar  to  Nusselt's 
original  film  condensation  film  theory  may  be  used  [Ref.  11] . 
Based  on  assumption  c,  a  static  force  balance  may  be  taken 
on  an  infinitesimal  fluid  element  in  the  x-direction  as  shown 
in  Figure  3.  This  force  balance  results  in  the  following 
equation: 

EFx  “  0  :  If  '  ll  =  0  Ceqn  2,1) 

where  x  =  shear  stress  (lbf/ft2) 
p  =  pressure  (lbf/ft2) 

x  =  co-ordinate  measuring  distance  along  surface  (ft) . 
y  =  co-ordinate  measuring  distance  normal  to  surface  (ft). 

3 .  Condensate  Momentum  Equation  (Y- Direction) 

In  a  similar  manner,  using  Figure  3,  a  force  balance 
in  the  y-direction  yields: 

EFv  =  0  :  ^2.  +  p  fu2r  =  0  (eqn  2.2) 

y  dy  1 

where  =  density  of  the  fluid  (lbm/ft3) 
oj  =  angular  velocity  (rad/hr) 
r  =  radius  (ft) 
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4 .  Fluid  Velocity 

Integrating  equation  (2.2)  between  the  limits  y  and 
S*  for  y  and  corresponding  limits  of  P  and  Py  for  pressure 
results  in  the  following  equation: 


P  =  Py  +  P£id*r  (S*-y) 


(eqn  2.3) 


where  Py  -  the  pressure  of  the  vapor  (lbf/ft2) ,  and 
6*  =  film  thickness  (ft) 

Differentiating  equation  (2.3)  with  respect  to  x  yields  the 
following  expression  for  dP/dx: 


dp 

dPy  + 

p£(jjzr  d6* 

dx 

dx 

dx 

(eqn  2.4) 


Applying  assumption  (f) ,  (Py  is  constant,  therefore,  dPv/dx=0) 
and  substituting  equation  (2.4)  into  equation  (2.1)  yields: 


3t 

dy 


pfoj2r 


dS* 

dx 


(eqn  2.5) 


Integrating  equation  (2.5)  with  the  corresponding  limits  of 
integration  y  to  5*  and  t  to  0  results  in  the  following  ex¬ 
pression  for  shear  stress: 


2  d6  *  r  » * 

x  =  pfu)  r  -  [y  -  6*] 

X  A  y 


(eqn  2.6) 
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U  U 

T  =  y  w 


(eqn  2.7) 


where  y  =  fluid  dynamic  viscosity  (lbm/ft-hr) 
u  =  condensate  velocity  (ft/hr) 

Substituting  equation  (2.7)  into  equation  (2.6)  and  integrating 
with  the  corresponding  limits  of  integration  0  to  y  and  0  to  u 


yields : 


pf“  r  d«*  f  y2 


S  IV  - 


(eqn  2.8) 


The  average  velocity  of  the  condensate  may  be  found  in  the 


following  manner: 


f  if 

u  =  /  ud y  =  -W 

6  O  K  J  r 


6  pfu)2r  d6*  p  2 


6  o  y  dx 


-  y<$*  dy  (eqn  2.9) 


d6Z 

u  dx  1  2  J 


5.  Continuity  Equation 


(eqn  2.10) 


The  continuity  equation  for  mass  flow  requires  that: 


m  =  P  £  u  A 


(eqn  2.11) 


where  m  =  condensate  mass  flow  rate  (lbm/hr) 

A  ■  cross  sectional  area  of  the  fluid  (ft2) 


This  can  also  be  written  as 


/*<5 

m  =  /  pfu  2irrdy 
Jo 


(eqn  2.12) 


Substituting  equation  (2.10)  into  equation  (2.12)  and  inte¬ 
grating  yields: 


2irpf  w2r2  d<5*  6* 

"  =  ' -  y -  dT  T 


m 


(eqn  2.13) 


Differentiating  this  equation  with  respect  to  x  yields: 


dm  _  _  2irpf  o)2r2  d  ,  dd  *  <$* 3 

dx  P  dx  dx  3 


(eqn  2.14) 


6 .  Energy  Equation 

.Having  applied  assumption  (e) ,  if  the  film  surface 
temperature  is  at  the  saturation  temperature  (Tsat)  of  the 
vapor  and  if  the  wall  of  the  axial  increment  is  at  a  given 
constant  temperature  (Tw) ,  then  the  heat  transfer  by  conduc¬ 
tion  of  a  fluid  element  of  surface  area  dA  is: 


dq  =  kf  (Tsat  ~  Tw)  dA 
6* 


(eqn  2.15) 


where  dq  =  differential  heat  transfer  rate  (Btu/hr) 
dA  =  2irrdx(ft2) 

k£  =  thermal  conductivity  of  the  condensate  film 


(Btu/hr- ft- F) 


T  0+.  =  saturation  temperature  (degrees  F) 

Sal 

T  =  inside  condenser  wall  temperature  (degrees  F) 

Considering  the  change  of  phase  and  defining  E£g  as 
the  average  enthalpy  change  of  the  vapor  in  condensing  to  a 
liquid  and  subcooling  to  the  average  liquid  temperature  of 
the  film,  then  dq  is  also  defined  by: 


dq  =  hfg  dm 


(eqn  2.16) 


where  h£g  =  laten  heat  of  vaporization  (Btu/lbm) 


Cp  =  specific  heat  (Btu/lbm  R) 


AT  =  (T 


sat 


-  T  ) 


vr 


hf„  =  hf„  +  0.35  •  c 
fg  fg  P 


AT 


Rearranging  equation  (2.16)  and  substituting  this  equation  into 
equation  (2.15)  yields: 


=  kf(Tsat  ~  Tw)  ^TTr 
dx  7* 


(eqn  2.17) 


Finally  coupling  the  energy  and  continuity  equations  result  in 
the  following  differential  equation: 


s*  4.  ,*’]  -■  3kf(Tsat  -  Tw)p 


dx  dx 


pf2ai2r  E£g 


(eqn  2.18) 


Equation  (2.18)  can  be  solved  using  the  Finite  Element 


Method  to  provide  the  film  thickness  profile  along  the  axial 


length  of  a  cylindrical  condenser.  Appendix  A  provides  a 
detailed  description  of  this  solution.  Once  the  film  pro¬ 
file  is  known,  a  steady  state  two-dimensional  heat  conduc¬ 
tion  analysis  can  be  performed. 

7 .  Determination  of  Heat  Transfer  Rate 

Assume  that  the  cylindrical  condenser  section  of  the 
rotating  heat  pipe  is  divided  axially  into  a  number  of  incre¬ 
ments.  Then  for  any  axial  increment  of  a  cylindrical  conden¬ 
ser,  the  differential  heat  flux  can  be  determined  by  the 
following  expression. 


dq 11  = 


(Tsat  *  Too) 

■  ■■  ■  ■— 

6  ^  thick  1 


(eqn  2.19) 


where  Too 

thick 

kw 


h 


ext 


ambient  temperature  (degrees  F) 
thickness  of  the  condenser  wall  (ft) 
thermal  conductivity  of  the  wall  material 
(Btu/hr- ft- F) 

external  heat  transfer  coefficient 
(Btu/hr- ft2 -F) 


Note  the  three  terms  in  the  denominator  are  the  thermal  re¬ 
sistances  of  the  film,  wall  and  external  convection  respectively. 

The  differential  heat  transfer  rate  for  any  increment 
can  be  found  by  the  following  relationships: 


dq  =  dq11  •  2irrdx 


(eqn  2.20) 
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or 


dq  =  -1  Iwl  1*. .  (eqn  2.21) 

<5  A  thick  J  1 

W  ’w" 

Equation  (2.21)  represents  the  total  heat  transfer  rate  for 
an  incremental  section  of  width  dx.  To  find  the  total  heat 
transfer  rate  for  the  entire  cylindrical  condenser,  the  in¬ 
cremental  heat  rates  must  be  summed  over  the  entire  length 
of  the  condenser.  Therefore: 

NDIV 

Qtotal  "  <*<1  te1n  2'22) 

where  NDIV  =  total  number  of  axial  increments. 

C.  THEORY  FOR  A  CYLINDRICAL  AXIALLY  FINNED  CONDENSER 
1 .  Assumptions 

Referring  to  Figure  4,  it  is  obvious  that  the  analysis 
of  a  cylindrical  internally  finned  condenser  is  more  compli¬ 
cated  due  to  the  mass  flow  from  the  fins  into  the  trough 
region  between  the  fins.  For  this  reason,  in  addition  to 
the  simplifying  assumptions  made  for  the  smooth  condenser 
which  are  listed  in  the  previous  section,  the  following 
assumptions  must  also  be  made: 

a)  Referring  to  Figure  5,  the  mass  flow  along  the  fin  sur¬ 
face  does  not  flow  axially  in  the  x-direction,  but  only 
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along  the  surface  of  the  fin  in  the  z-direction  into  the 
trough.  Thus,  mass  flow  in  the  axial  direction  is  only- 
permitted  in  the  trough  region  between  the  fins.  This 
is  a  reasonable  assumption  in  that  the  film  thickness 
along  the  fin  surface  is  very  small  in  relation  to  the 
film  thickness  in  the  trough.  This  being  the  case,  the 
hydrostatic  force  in  the  x-direction  on  a  fin  fluid  ele¬ 
ment  will  be  much  less  than  the  centrifugal  force  com¬ 
ponent  in  the  z-direction  on  that  same  fluid  element 
forcing  that  fluid  element  into  the  trough. 

b)  Just  as  in  the  axial  direction,  there  is  no  pressure 
change  along  the  surface  of  the  fin  in  the  z-direction. 

c)  It  will  be  assumed  that  the  temperature  along  the  con¬ 
vective  surface  of  the  fin  is  at  a  constant  value  CTavg) 
This  average  fin  surface  temperature  is  the  arithmetic 
average  of  the  fin  tip  temperature  and  the  fin  base  sur¬ 
face  temperature  where  the  fin  intersects  with  the  wall 
of  the  condenser.  This  is  a  valid  assumption  if  the  fin 
section  is  divided  into  a  sufficient  number  of  finite 
elements.  Purnomo’s  [Ref.  5]  results  indicate  a  less 
than  one  degree  variation,  even  for  very  large  fin  half 
angles.  This  variation  in  temperature  will  have  an  in¬ 
significant  effect  on  film  thickness  along  the  surface 
of  the  fin  and  can  be  neglected  by  using  an  average 


value . 


Mass  Flow  in  the  X-Direction 


As  a  result  of  assumption  (a) ,  the  resulting  momentum 
equation  in  both  the  x  and  y  directions  as  well  as  the  equa¬ 
tions  for  velocity  and  mean  velocity  are  identical  to  those 
developed  in  section  B  for  the  smooth  condenser  and  will  not 
be  redeveloped  here.  Looking  now  at  mass  flow  in  the  x-direc- 
tion  which  is  limited  to  flow  only  in  the  trough,  the  mass 
flow  rate  is  given  by  the  following  expression: 


"'total 


d6* 

dx 


<5*2(e5*  +  5* 2  tana) 


(eqn  2.23) 


where  a  =  fin  half  angle  (radians) 
e  =  width  of  the  trough  (ft) 

Note  that  the  quantity  in  parentheses  is  the  cross  sectional 
area  of  the  film  condensate  in  the  trough  (See  Figure  4) . 
Taking  the  derivative  of  equation  (2.23)  with  respect  to  x 
yields  the  rate  of  change  of  mass  flow  in  the  trough  for  a 
given  axial  increment. 


dl"total 

dx 


Pfu)2r 

3y 


—  iff  (eo*3+  5*  4tana) ] 
dx  dx 


(eqn  2.24) 


Equation  (2.24)  represents  the  rate  of  change  of  the 
total  mass  flow  rate  with  respect  to  x  in  the  x-direction. 
This  equation  must  be  coupled  with  the  energy  equations  for 
the  fin  and  trough  to  develop  a  representation  of  the  film 
profile  in  the  trough. 


3.  Mass  Flow  in  the  Z-Direction 


Examining  an  infinitesimal  fluid  element  on  the  sur¬ 
face  of  a  fin  for  any  axial  increment  of  width  Ax,  as  shown 
in  Figure  6,  the  momentum  equation  in  the  z-direction  becomes: 


3P 

3z 


P£(u2rcosa 


(eqn  2.25) 


where  tz  =  shear  stress  in  the  z-direction  (lbf/ft2) 

z  =  co-ordinate  measuring  distance  along  the  surface 
of  the  fin  (ft) 

Neglecting  dP/dz  based  on  assumption  (b) >  and  integrating 
equation  (2.25)  from  t  to  0  and  y  to  6  yields: 

t2  -  u  -  pfw2rcosa(S-y)  (eqn  2.26) 


where  w  =  fiuid  velocity  in  the  z-direction  (ft/hr) 

6  =  fin  film  thickness  along  the  surface  of  the  fin  (ft) 
Note,  6  ,  the  fin  film  thickness  should  not  be  confused  with 
<5*,  the  film  thickness  in  the  trough.  Integrating  equation 
(2.26)  from  0  to  w  and  0  to  y  provides  the  following  expression 
for  fluid  velocity: 

pf62rcosa  2 

W  -  - -  (S(z)y  -  \  )  (eqn  2.27) 


This  relationship  may  be  used  to  find  the  average  fin  fluid 
velocity  w: 


p^to2rcosa  <S(z)2 
3m 


(eqn  2.28) 


The  mass  flow  rate  in  the  z-direction  along  the  surface  of  the 
fin  for  a  given  axial  increment  is  given  by: 

rii£^n  =  P  £wdA  (eqn  (2.29) 

where  dA  =  the  cross  sectional  area  of  the  fluid  flowing  along 
the  fin  surface  (ft2) 

Substituting  equation  (2.28)  into  equation  (2.29)  yields: 
pf2oj2rcosa  63(z)dx 

mfin  =  -  (eqn  2.30) 

3y 

This  equation  is  identical  to  equation  (15)  of  Schafer's  [Ref. 

3]  analysis  if  the  condenser  cone  half  angle  (0)  is  set  equal 
to  0  which  is  the  case  for  a  cylindrical  condenser. 

4 .  Energy  Equation  for  the  Trough  Condensate 

An  energy  balance  on  an  infinitesimal  fluid  element  in 
the  trough  of  an  axial  increment  of  width  dx  with  surface  area 
e  •  dx  yields  the  following  expression  for  heat  transfer  by 
conduction: 


-  TJedx 


+■  rnn  rrVi 


-  kf<Tsat 


(eqn  2.31) 


Note  also,  that  the  trough  heat  transfer  rate  is  given  by: 


^trough  "  hfg  dAtrough  (ec>n  2' 

Combining  equations  (2.31)  and  (2.32)  and  dividing  by  dx  re¬ 
sults  in  the  following: 


^trough  _  ^f^Tsat  "  Tw^c 
dx  Efg  6* 


(eqn  2.33) 


Equation  (2.33)  is  an  expression  for  incremental  change  in  mass 
flow  rate  with  respect  to  x  due  to  condensation  in  the  trough 
region. 

5 .  Energy  Equation  for  the  Fin  Condensate 

An  energy  balance  on  a  differential  element  of  surface 
area  dx*dz  yields  the  following  relationship  for  differential 
heat  into  the  fin: 


d«£in 


hfg  d*f  in 


kf<Tsat 


-  Tfin(z))dxdz 
5(z) 


(eqn  2.34) 


where  T^n  (z)  =  fin  surface  temperature  at  some  position  z 

along  the  surface  of  the  fin  (degrees  F) 

Since  the  fin  condensate  mass  is  assumed  to  flow  only  in  the 
z-direction,  equation  (2.30)  is  differentiated  with  respect 
to  z  and  substituted  into  equation  (2.34).  After  substitution 
and  rearrangement,  the  following  equation  results: 
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5  (z)  3d<5  (z)  = 


VTsat  -  W2»da 

P£2oi2r  h^g  cosa 


(eqn  2.35) 


Applying  assumption  (c)  i.e.,  Tfin(z)  equals  Tavg  for 
all  z  and  integrating  equation  (2.35)  from  0  to  6  and  0  to  z 
yields  the  following  relationship  for  fin  film  thickness  <5  (z)  : 

i  A 


6(z)  = 


4  MT 


sat 


-  T  )  y  z 
avgJ 


P£2ui2r  h^g  cosa 


(eqn  2.36) 


where  Tavg  =  average  fin  surface  temperature  (degrees  F) . 
Substituting  equation  (2.36)  into  equation  (2.30)  and  solving 
for  rate  of  change  of  mass  flow  rate  of  the  fin  with  respect 
to  x  for  an  increment  of  width  dx  yields: 


dlilfin 

dx 


2Pf 2ai2rcosa  f4  kf(Tsat  -  Tavg)yz* 


3y 


■u2r  h 


fg 


cosa 


(eqn  2.37) 


where  z*  =  z-6*/cos(a).  Note,  z*  is  the  distance  along  the 
surface  of  the  fin  from  the  fin  tip  to  the  trough  film  thick¬ 
ness  (6*).  Note  also,  that  the  right  hand  side  of  equation 
(2.37)  is  multiplied  by  two;  this  accounts  for  mass  flow  from 
the  fins  on  both  sides  of  the  trough. 

6 .  Continuity  Equation 

For  any  axial  increment  of  length  dx,  continuity  dic¬ 
tates  that: 
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total  =  fin  +  umt rough 
3x  dx  3x 


(eqn  2.38) 


Substituting  equations  (2.24), 
(2.38)  and  rearranging  yields: 

5*  4  [aT  Cs6*3+  5*4tana^  - 


(2.33)  and  (2.37)  into  equation 


3  kffTsat  -  T«>^ 
pf2“2r  Efa 


2  5  *cosa 


4  kf<Tsat 


-  T  )yz* 
avg-'  M 


3/4 


>£  wzr  h£g  cosa 


(eqn  2.39) 


Equation  (2.39)  can  be  solved  using  the  Finite  Element  Method 
formulation  provided  in  Appendix  A.  The  solution  of  this  equa¬ 
tion  provides  the  film  thickness  profile  along  the  axial  length 
of  a  cylindrical  finned  condenser. 

7 .  Determination  of  the  Heat  Transfer  Rate 

Once  the  film  profile  has  been  determined  within  the 
trough,  the  local  convective  heat  transfer  coefficient  can  be 
found  for  the  trough  using  the  following  relationship: 


h(x)trough  *  ^*Yx)  Ceqn  2'40:i 

In  a  similar  manner,  the  local  heat  transfer  coefficient  along 
the  surface  of  the  fin  can  be  found  by: 


h(z) 


fin 


5 


(eqn  2.41) 


The  differential  heat  transfer  rate  for  any  fin  section,  as 
shown  in  Figure  4,  of  axial  incremental  length  dx  is: 
a)  for  the  trough: 


dq 


trough 


^sat  -  Vzdx 
trough 


(eqn  2.42) 


where  e*dx  is  the  surface  area  of  the  trough,  and 
b)  for  the  fin  surface: 


dq 


fin 


CTsat  -  T 


h(z) 


avg)dxdz 

fin 


(eqn  2.43) 


where  zq  is  the  surface  length  of  the  fin. 

The  total  differential  heat  transfer  rate  per  axial  increment 
is  found  by  summing  equation  (2.42)  and  (2.43)  for  the  total 
number  of  fins.  That  is: 


NFIN 

d«total  =  l  Cdl>fin  *  ^trough3  (e£In  2’44) 


where  NFIN  is  the  total  number  of  axial  fins. 

In  a  similar  manner,  the  total  heat  transfer  rate  for  the  entire 
finned  condenser  can  be  found  by  the  following  relationship: 

NDIV 

^total  ^  ^total  (eqn  2.45) 


— ( 


where  NDIV  is  the  total  number  of  axial  increments. 


III.  COMPUTER  CODE  DESCRIPTION 


A.  GENERAL  DESCRIPTION  OF  CODE 

The  computer  code  consists  of  a  main  body  and  eight  sub¬ 
routines.  Basically,  the  code  which  is  provided  in  Appendix 
C  is  a  modification  of  Pumomo's  [Ref.  5]  code.  The  function 
of  each  subroutine  used  in  the  code  is  as  follows: 

a)  MCORRESM  established  the  correspondence  between  the 
local  and  global  nodal  points  used  in  the  Finite  Element 
Method  solution  for  the  two-dimensional  steady  state 
heat  conduction  problem.  In  so  doing,  "CORRES"  also 
numbers  all  elements  and  nodal  points  in  the  finite  ele¬ 
ment  model  and  assigns  local  nodal  points  to  each  of  the 
elements.  In  addition  MCORRESM  also  defines  major  ele¬ 
ment  numbers  used  in  other  subroutines  as  control 
parameters . 

b)  "COORD"  defines  the  x  and  y  coordinates  for  all  nodal 
points  in  the  finite  element  heat  conduction  problem 
model. 

c)  "DLSTAR"  determines  the  film  thickness  (6*)  on  the  sur¬ 
face  of  a  smooth  condenser  or  in  the  trough  in  a  finned 
condenser . 

d)  MHTCOEF"  determines  the  heat  transfer  coefficient  for 
all  convective  surface  elements. 


e)  "FORMAF"  formulates  the  Finite  Element  Method  equations 
for  the  two-dimensional  steady  state  heat  conduction 
problem. 

f)  "BANDEC"  is  an  equation  solver  for  a  symmetric  matrix 
which  has  been  transformed  into  banded  form.  "BANDEC" 
will  return  the  solution  to  the  two-dimensional  heat 
conduction  problem. 

g)  "HTCALC"  determines  the  elemental,  incremental  and  total 
heat  transfer  rates. 

h)  "DELCRV"  determines  the  condensate  film  profile  in  a  cyl¬ 
indrical  condenser. 

Two  additional  Naval  Postgraduate  School  computer  library 
routines  are  also  used  in  the  code: 

a)  "DPOLRT"  is  a  nonIMSL  double  precision  library  routine 
that  determines  the  roots  of  a  real  polynomial.  This 
routine  is  called  by  "DLSTAR"  to  determine  the  film 
thickness  for  the  succeeding  increment  in  the  analysis 
of  a  tapered  condenser. 

b)  "LEQT2F"  is  an  IMSL  double  precision  libary  routine  that 
solves  a  set  of  simultaneous  linear  equations.  This 
routine  is  called  by  "DELCRV"  to  solve  the  Finite  Element 
Method  equations  for  the  cylindrical  condenser  film  pro¬ 
file  problem.  The  resulting  film  profile  is  then  used 

in  the  heat  conduction  analysis. 

In  order  to  use  the  computer  code  to  analyze  heat  transfer 
in  a  rotating  heat  pipe,  nine  data  cards  are  required.  A  user's 
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guide  describing  these  data  cards  and  required  input  is  pro¬ 
vided  in  Appendix  B.  The  input  data,  describing  the  geometric 
configuration  of  the  rotating  heat  pipe  as  well  as  the  operating 
parameters  determines  which  solution  technique  is  utilized  in  the 
analysis.  The  solution  technique  for  each  of  the  four  con¬ 
denser  geometries,  i.e.,  tapered-smooth ,  tapered-axially  finned, 
cylindrical-smooth  and  cylindrical-axially  finned  is  different. 

In  all  cases  however,  the  Finite  Element  Method  is  used  to 
solve  the  two-dimensional  steady  state  heat  conduction  problem. 
This  solution  is  the  one  developed  by  Purnomo  [Ref.  5]  and 
has  not  been  modified.  Details  of  the  development  of  this 
solution  are  described  in  detail  in  Purnomo' s  thesis  [Ref.  5] 
and  will  not  be  repeated  here.  This  being  the  case,  each  of 
the  four  solution  techniques  will  now  be  discussed  in  detail. 

B.  INTERNALLY  FINNED  TAPERED  CONDENSER  SOLUTION 

The  complete  development  of  this  solution  technique  is 
provided  in  Purnomo 's  [Ref.  5]  thesis  and  will  not  be  rede¬ 
veloped.  When  an  equation  is  required  for  clarity,  the  equa¬ 
tion  in  final  form  will  be  provided.  Where  there  is  a 
modification  to  Purnomo 's  [Ref.  5]  code,  this  modification 
will  be  noted. 

The  condenser  of  the  internally  finned- tapered  condenser 
is  divided  into  NDIV  axial  increments.  These  axial  incre¬ 
ments  are  then  subdivided  circumferentially  into  ZFIN  number 
of  subincrements  where  ZFIN  is  the  total  number  of  fins. 
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These  subincrements  are  then  divided  in  half  to  form  the 
basic  symmetric  unit,  one-half  of  a  fin-trough  section  as 
shown  in  Figure  1.  This  unit  section  is  then  divided  into 
a  number  of  linear  triangular  elements.  The  number  of  ele¬ 
ments  depend  on  the  input  parameters.  The  only  limitation 
is  that  the  number  of  system  nodal  points  must  hot  exceed 
100;  otherwise,  certain  variables,  e.g.,  x  and  y,  would 
exceed  allotted  storage  values.  Figure  7  shows  a  unit 
section  subdivided  into  25  elements.  After  this  unit  is 
subdivided,  each  nodal  point  is  assigned  an  x  and  y  coordi¬ 
nated  based  on  the  geometric  input  parameters. 

To  start  the  iteration,  two  initial  values  are  required: 
1)  an  initial  temperature  for  the  nodal  points  along  the 
internal  convective  boundary,  and  2)  an  initial  trough  film 
thickness  (6*).  The  initial  temperature  is  provided  as  an 
input  parameter  and  the  initial  trough  film  thickness  at  the 
first  increment  is  provided  by  a  relationship  taken  from  an 
analysis  by  Sparrow  and  Gregg  for  condensation  on  a  rotating 
disk  [Ref.  12 ] . 

Once  these  values  are  known,  the  heat  transfer  coeffi¬ 
cient  for  the  internal  convective  surface  elements  are  found 
using  the  following  relationships: 


i/ * 


h(z) 


=  ^f 

fin  <ff(  z) 


■[ 


kfPf2wz  (r+xsin0  )h£gcos0cosa 
4u(-AA.z373-BB-z2/2+CTsat-T1)z 


(eqn  3.1) 
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Figure  7.  Axially  Finned  Condenser  Symmetric  Section  Sub¬ 
divided  into  25  Linear  Triangular  Finite  Elements 


where  x  =  distance  from  condenser  end  (x=0.0)  to  midpoint 
of  increment  (ft) 

T-^  =  fin  apex  temperature  (degrees  F) 

AA,  BB  =  constants  in  the  parabolic  temperature  determined 
by  a  Langrangian  fit. 

For  the  trough  surface  elements : 

htx)trough  *  pfx)  teqn  3'2) 

These  heat  transfer  coefficients,  along  with  the  thermal 
conductivity  and  x  and  y  coordinates  of  the  nodal  points  are 
used  to  form  the  Finite  Element  Method  equations  for  the  two- 
dimensional  heat  conduction  problem.  The  equations  are  then 
solved  to  yield  a  temperature  distribution  in  the  symmetric 
section. 

The  above  iteration  is  repeated,  where,  now  the  solution 
temperature  distribution  from  the  previous  iteration  is  used 
to  calculate  the  heat  transfer  coefficients  along  the  con¬ 
vective  surface  fin  elements  as  well  as  a  new  5*  using  Sparrow's 
and  Gregg's  relationship  [Ref.  12).  This  new  5*,  in  turn,  is 
used  to  determine  the  heat  transfer  coefficients  of  the  trough 
elements . 

Again,  the  Finite  Element  Method  subroutines  will  yield 
a  temperature  distribution  for  the  symmetric  section.  At 
this  point,  the  nodal  point  temperatures  are  checked  for  con¬ 
vergence  using  the  following  relationship: 


Max 


<  CRIT  i  =  1,2, . . . NSNP 


(eqn  3.3) 


T-  .  -  T.  .  , 
T-  • 


where  NSNP  =  number  of  system  nodal  points 
j  *  present  iteration 
j-1  =  previous  iteration 
CRIT  =  convergence  criterion 

If  this  convergence  test  is  successful  for  all  nodal  points 
the  increment  is  considered  solved.  If  any  one  nodal  point 
fails,  the  iterative  process  is  repeated  until  convergence  is 
met.  The  convergence  test  in  equation  C 3 . 3)  is  different  than 
the  one  used  by  Purnomo  [Ref.  5]  in  his  thesis.  Purnomo  com¬ 
pared  incremental  heat  transfer  rates  per  unit  of  condenser 
length,  Q^,  rather  than  temperature  as  is  done  in  the  modified 
code . 

If  convergence  is  met,  the  heat  transfer  rate  is  determined 
From  this  heat  transfer  rate,  the  incremental  mass  flow  rate  is 
determined  by  the  following  equation: 


"‘total 


2Qt  Ax 

~1^T 


(eqn  3.4) 


where  =  heat  transfer  rate  per  unit  length  (Btu/hr-ft) 

Ax  =  incremental  width  (ft) 

Using  this  value  of  incremental  mass  flow  rate  determined 
by  equation  (3.4),  the  following  equation  is  used  to  calculate 
the  subsequent  interval's  trough  condensate  film  thickness  (<$*) 
with  a  polynomial  rootfinder  subroutine: 
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®total 


Pf2w2  (r+xsinj3)  6  *2  (x)sin0(<S*(x)  e  +  6*2  (x)  tana) 


3y 


(eqn  3.5) 


where  e  =  trough  width  (ft) 

This  resulting  value  of  6*(x)  is  then  defined  as  the  trough 
film  thickness  for  the  next  increment.  In  addition,  the  solu¬ 
tion  temperature  distribution  from  the  previous  iteration  is 
used  as  the  starting  temperature  distribution  for  the  next 
increment . 

This  iterative  process  at  each  increment  is  repeated  until 
convergence  is  met,  and,  is  continued  at  each  increment  until 
the  entire  length  of  the  condenser  has  been  transversed.  In¬ 
cremental  heat  rates  are  then  summed  to  yield  the  total  heat 
transfer  rate.  That  is: 

NDIV 

Qtotal  =  2*ZFIN*  l  Q^Ax  •  (eqn  3.6) 

where  ZFIN  is  the  total  number  of  axial  fins. 

Once  the  total  heat  transfer  rate  has  been  determined,  the 
problem  is  solved  and  pertinent  data  is  provided  as  output. 

C.  SMOOTH  TAPERED  CONDENSER  SOLUTION 

The  heat  pipe  condenser  is  divided  into  NDIV  number  of 
axial  increments  as  in  the  finned- tapered  condenser  solution. 
These  axial  increments  are  then  subdivided  into  360  segments 
of  equal  length;  these  segments  are  the  basic  symmetric  unit. 
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This  unit  section,  is  divided  into  a  number  of  linear  tri¬ 
angular  elements  with  the  same  limitation  as  before;  the 
number  of  system  nodal  points  must  not  exceed  100.  The 
system  nodal  points  are  then  assigned  x  and  y  coordinates 
based  on  the  input  geometric  parameters. 

To  start  the  iterative  process,  as  in  the  finned- tapered 
case,  the  initial  value  of  temperature  which  is  an  input 
parameter  is  used  to  solve  for  the ‘ initial  value  of  fin 
thickness  (6*)  based  on  the  Sparrow  and  Gregg  analysis  [Ref. 
12]  . 

Once  this  initial  value  of  6*  is  known,  the  heat  transfer 
coefficients  for  the  internal  convective  elements  can  be 
determined  using  equation  (3.2).  These  heat  transfer  coef¬ 
ficients  are  used  in  the  Finite  Element  Method  equations. 

The  equations  are  solved  yielding  a  temperature  distribution. 
The  iteration  is  repeated  until  convergence  is  met,  just  as 
in  the  finned- tapered  case. 

When  convergence  is  met,  that  is  equation  (3.3)  has  been 
satisfied,  a  new  film  thickness  6*(x)  can  be  found  by  one  of 
the  following  equations  for  6*(x): 

Shx[V]  '  T  Drx  [V]  ‘  T  Revx  [1T]  '  1  =  0  (ec*n  3-7) 

or 


(eqn  3.8) 
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S*(x)  = 


*V<Tsat-V 

2p£2w2rsin20h£g 


[1  -  {- 


(r+xsin0) 


73 
-}  ] 


1  /  4 


(eqn  3.9) 


where  Sh  = 


pf2  (a>2r-g)sin0E^gX; 
4^MTsat  -  V 


Dr  = 


pftvhfgx2cos0 

"MTsat  -  V 


Rev  = 


P£V  cos0 


g  =  acceleration  due  to  gravity  (ft/hr2) 
v  =  vapor  velocity  (ft/hr) 

Tw  =  local  wall  temperature  (degrees  FO 

=  shear  stress  vapor-liquid  interface  (lbf/ft2) 
Equation  (3.7)  defines  the  film  thickness  distribution  for 
a  smooth  tapered  rotating  heat  pipe  derived  by  Daniels  and 
Al-Jumaily  [Ref.  13].  This  equation  takes  into  account  the 
drag  effects  of  counter- flowing  vapor.  Equation  (3.8)  is  a 
modification  of  Equation  (3.7)  neglecting  the  drag  losses. 
Equation  (3.9)  was  developed  by  Ballback  [Ref.  1].  This  equa¬ 
tion  also  neglects  drag. 

Depending  on  a  particular  control  parameter  which  is  part 
of  the  input  data  (See  Appendix  B)  one  of  these  equations  is 
used  to  solve  for  the  film  thickness  (6*(x))  for  the  next 
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increment.  In  addition,  the  solution  temperature  distribution 
from  the  previous  iteration  is  used  as  the  starting  temperature 
distribution  for  the  next  increment. 

This  iterative  process  at  each  increment  is  repeated  until 
convergence  is  met,  and  is  continued  at  each  increment  until 
the  entire  length  of  the  condenser  has  been  transversed,  just 
as  in  the  f inned-tapered  case.  Total  heat  transfer  rates  are 
then  determined  by  summing  the  incremental  heat  transfer  rates 
for  the  entire  length  of  the  condenser  by  the  following 
relationship : 

NDIV 

Qtotal  *  36°*  .Z  Qj.Ax  Oqn  3*1J) 

D.  SMOOTH  CYLINDRICAL  CONDENSER  SOLUTION 

As  in  the  smooth- tapered  case,  the  condenser  is  first 
divided  axially,  then  it  is  divided  circumferentially  into 
360  segments  of  equal  length.  These  segments  are  the  basic 
symmetric  unit  section  to  be  considered.  Again  the  symmetric 
unit  is  subdivided  into  linear  triangular  elements  and  x  and 
y  coordinates  are  assigned  to  the  system  nodal  points. 

To  begin  the  iteration,  an  initial  temperature  estimate 
which  is  an  input  parameter  is  assigned  to  the  convective  sur¬ 
face  nodal  points.  Using  this  initial  temperature  estimate, 

the  maximum  film  thickness,  6*  „  which  is  located  at  x  =  0, 

max 

is  determined.  This  maximum  film  thickness  value  is  then 


used  as  one  of  the  boundary  conditions  in  the  solution  of 
equation  (2.18)  using  the  Finite  Element  Method.  Equation 
(2.18)  is  repeated  here  for  reference. 


5*  Sx  4r  «*’)  -  _  V  (eqn  3.11) 

Pf  hfg 

The  finite  element  solution  of  Equation  (3.11)  provides  the 
film  thickness  profile  (5*(x))  at  the  midpoint  of  each  incre¬ 
ment  along  the  length  of  the  condenser,  but  only  the  first 
increment  value  is  applied  at  this  point  of  the  analysis. 

Once  <5*(x)  is  known  at  the  first  increment  the  heat  transfer 
coefficients  for  the  internal  convective  surface  elements  can 
be  determined  using  equation  (3.2).  The  steady  state  heat  con¬ 
duction  problem  is  then  solved  and  a  temperature  distribution 
for  the  unit  section  results.  This  process  is  now  repeated. 

A  new  maximum  film  thickness,  film  profile  and  thus  5*(x)  at 
the  first  increment  is  found  based  on  the  solution  temperature 
distribution  from  the  first  iteration.  With  this  new  value  for 
(x) ,  the  heat  conduction  problem  is  again  solved  for  a  new 
temperature  distribution.  This  iterative  process,  involving 
solution  of  the  film  profile  with  each  iteration  is  continued 
until  temperature  convergence  is  met  at  the  first  increment. 
When  convergence  is  met,  the  film  profile  determined  on  the 
iteration  in  which  convergence  was  met  is  used  to  provide  the 
values  of  5*(x)  for  the  remaining  increments  along  the  length 
of  the  condenser. 


Using  this  predetermined  value  of  6*(x),  the  iterative  pro¬ 
cess  is  continued  at  each  increment  until  temperature  conver¬ 
gence  is  reached.  As  convergence  is  reached  at  each  increment, 
the  internal  wall  temperature,  which  is  the  same  for  all  inter¬ 
nal  wall  nodal  points  of  the  section,  is  stored  for  future 
application. 

When  convergence  is  reached  at  the  final  increment,  equation 
(3.11)  is  once  again  solved  for  the  film  profile.  It  should  be 
noted,  however,  that  the  right  hand  side  of  equation  (3.11)  is 
temperature  dependent.  It  should  also  be  noted  that  the  temper¬ 
ature  varies  axially  along  the  length  of  the  condenser.  In 
order  to  account  for  this  temperature  dependence  and  temperature 
variation,  the  right  side  of  equation  (3.11)  is  now  determined 
for  each  increment  using  the  wall  temperatures  that  were  stored 
at  each  increment.  The  finite  element  solution  of  the  film  pro¬ 
file  will  now  account  for  the  temperature  variation  along  the 
length  of  the  condenser. 

This  final  film  profile  provides  the  value  of  S*(x)  for  each 
increment.  The  iterative  process  of  solving  for  heat  transfer 
coefficient,  temperature  distribution,  and  temperature  conver¬ 
gence  is  continued  for  each  increment  until  the  final  increment 
is  reached. 

Equation  (2.13)  provides  a  relationship  for  total  mass  flow 
rate  as  a  function  of  position.  The  total  mass  flow  rate  at 
the  overfall  into  the  evaporator  is  given  by: 


(eqn  3.12) 


m 


2irpf  2g)2r  .  d6*(L) 

p  dx 


3* 3 (L) 
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where  S*(L)  and  d5*(L)/dx  are  the  values  for  the  film  thickness 
and  rate  of  change  of  film  thickness  with  respect  to  x  at  the 
evaporator  end  of  the  condenser.  Another  relationship  for  mass 
flow  rate  based  on  ‘the  steady  state  heat  conduction  solution 
is  given  by: 


m  = 


NDIV 
360*  Z 
i=l 


"fg 


(eqn  3.13) 


If,  in  fact,  the  solution  of  equations  (3.12)  and  (3.13)  are 
equal,  then  the  mass  flow  rate  of  the  condensate  returning  to 
the  evaporator  is  equal  to  the  mass  flow  rate  of  the  vapor 
being  condensed  on  the  surface  of  the  condenser.  Or,  to  put 
it  another  way,  continuity  is  satisfied. 

It  should  be  noted  that  the  film  profile  maintains  the  same 
basic  shape,  that  is,  5*(x)  at  x=0  is  always  equal  to  5*ax  and 
decreases  to  a  specific  minimum  value  at  x=L.  This  being  the 
case,  if  the  maximum  film  thickness  is  varied,  the  film  thick- 
nexx  profile  will  vary  in  the  same  manner.  For  example,  if 
the  maximum  film  thickness  is  increased,  the  entire  profile 
will  also  increase.  This  will  result  in  an  increased  internal 
thermal  resistance  and  thus  lower  heat  transfer  rate.  As  a 
result  of  the  lower  heat  transfer  rate,  the  mass  flow  rate  as 
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determined  by  equation  (3.13)  will  be  less.  At  the  same  time, 
however,  the  greater  film  profile  will  result  in  a  greater 
value  of  6*(x)  at  the  overfall,  x  =  L.  Yet,  since  the  profile 
maintained  the  same  basic  shape,  the  derivative  at  the  overfall 
remains  relatively  constant.  Thus  the  mass  flow  rate  as  deter¬ 
mined  by  equation  (3.12)  will  increase.  A  decrease  in  the 
maximum  film  thickness  will  result  in  an  opposite  effect  to 
the  mass  flow  rates  determined  by  equations  (3.12)  and  (3.13). 

This  being  the  case,  the  mass  flow  rates,  as  determined 
by  equations  (3.12)  and  (3.13)  are  now  compared  to  determine 
if  the  film  profile  is  in  fact  the  solution  profile  to  the 
problem.  If  continuity  is  not  satisfied,  the  maximum  film  thick¬ 
ness  is  varied  and  the  entire  iterative  process  is  restarted. 

This  process  is  continued  until  the  film  profile  mass  flow  rate, 
equation  (3.12)  converges  towards  the  heat  transfer  mass  flow 
rate,  equation  (3.13).  When  the  absolute  difference  between 
these  mass  flow  rates  is  less  than  a  specific  value,  the  re¬ 
sulting  heat  transfer  rate  is  considered  the  solution  to  the 
problem. 

E.  FINNED  CYLINDRICAL  CONDENSED  SOLUTION 

As  in  the  finned- tapered  case,  the  condenser  is  divided 
axially  and  circumferentially  into  the  basic  symmetric  section 
as  shown  in  Figure  1.  This  unit  is  then  subdivided  into  linear 
triangular  elements  and  x  and  y  coordinates  are  assigned  to 
each  nodal  point. 
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To  begin  the  iterative  process,  as  before,  an  initial  temp¬ 
erature  estimate  is  assigned  to  each  nodal  point  along  the  in¬ 
ternal  convective  surface.  An  initial  trough  film  profile  is 
determined,  using  equation  (3.11)  and  this  initial  temperature 
estimate.  In  this  case,  however,  the  maximum  film  thickness 
(5*ax)  is  not  calculated  but  is  an  input  parameter. 

Once  5*(x)  is  known  at  the  first  increment,  the  internal 
heat  transfer  coefficients  are  determined,  using  equations  (3.1) 
and  (3.2).  These  values  are  then  used  in  the  Finite  Element 
Method  solution  of  the  steady  state  heat  conduction  problem. 

A  temperature  distribution  is  determined  and  the  iteration  is 
repeated  until  temperature  convergence  is  met.  Note  that  a 
new  film  profile  is  determined  for  each  iteration. 

As  in  the  cylindrical-smooth  condenser  case,  once  conver¬ 
gence  is  met  at  the  first  increment,  the  film  profile  that  was 
determined  for  the  iteration  prior  to  convergence  at  the  first 
increment  is  then  used  to  provide  the  film  thickness  5*(x)  for 
the  remaining  increments. 

At  the  final  increment,  a  new  film  profile  for  a  finned 
cylindrical  condenser  is  then  determined  by  solving  the  follow¬ 
ing  equation  developed  in  Chapter  II: 


5*  (5*3e+  5 *4 tana)  ]  =  -  3kf(?*at,  I  Iw)ue 
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(eqn  3.14) 
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A  description  of  the  solution  of  this  differential  equation 
using  the  Finite  Element  Method  is  provided  in  Appenxix  A. 


The  iterative  process  of  finding  the  solution  temperature 
distribution  for  all  increments  is  then  repeated  until  the 
length  of  the  condenser  has  been  transversed. 

At  this  point,  the  total  mass  flow  rate  of  the  condenser 
returning  to  the  evaporator  is  determined  by  the  following 
relationship: 

Pf2W2r  Jr* 

m  =  -  ^ (e5*  3+6*4tana)  *  ZFIN  (eqn  3.15) 


where  ZFIN  is  the  number  of  axial  fins. 

This  mass  flow  rate  is  compared  to  the  mass  flow  rate  given 
by  the  following  equation: 


m  = 


ZFIN*2 


NDIV 
*  l 
i=l 


(eqn  3.16) 


Just  as  in  the  smooth-cylindrical  condenser  case,  if  the 
absolute  difference  between  the  two  mass  flow  rates  is  less 
than  a  mass  flow  convergence  criterion,  the  problem  is  con¬ 
sidered  solved.  If  not,  <$*  „  is  varied  and  the  entire  itera- 

max 

tive  process  is  started  again.  As  in  the  smooth-cylindrical 
condenser,  varying  s*ax  will  have  the  same  effect  on  the  film 
profile  and  the  heat  transfer  rate.  Since  the  temperature  dis¬ 
tribution  along  the  condenser  has  been  solved  once  and  closely 
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approximates  the  final  solution  to  the  problem,  on  successive 
iterations,  equation  (3.14)  may  be  used  to  solve  for  the  film 
profile  rather  than  equation  (3.11).  Equation  (3.11)  was  used 
on  the  first  iteration  because  the  finite  element  solution  con¬ 
verges  more  quickly  than  equation  (3.14)  when  an  estimated 
temperature  is  used. 

A  word  of  caution  is  required.  The  solution  of  equation 

(3.14)  is  highly  sensitive  to  the  value  of  6*ax-  If  the  initial 

value  of  6*  is  inconsistent  with  the  actual  solution,  e.g. 
max  ’ 

too  small,  the  Finite  Element  Method  solution  of  the  film  pro¬ 
file  will  not  converge.  If  this  is  in  fact  the  case,  the  prob¬ 
lem  will  be  automatically  terminated.  A  new  value  of  6*  should 

max 

then  be  chosen  and  the  problem  restarted. 

One  additional  topic  which  should  be  addressed  is  the  rec¬ 
tangular  fin  solution  process.  The  rectangular  fin  profile  is 
a  slight  modification  of  the  f inned- tapered  or  finned-cylindrical 
condenser  solution.  The  only  variation  is  that  the  top  elements 
of  the  rectangular  fin  are  assigned  heat  transfer  coefficients 
of  0.0.  Thus,  the  tip  of  the  rectangular  fin  is  considered 
adiabatic  and  no  heat  is  transferred  through  this  top  face. 

Other  than  this  modification,  the  solution  techniques  are  the 
same  as  for  the  finned  cases  addressed  above. 


IV.  RESULTS  AND  DISCUSSION 


Prior  to  the  heat  transfer  analysis  of  the  various  con¬ 
denser  configurations,  it  was  necessary  to  verify  the  finite 
element  solution  of  the  film  profile.  The  development  of  this 
solution  is  discussed  in  Appendix  A. 

Leppert  and  Nimmo  [Refs.  8  and  9]  had  developed  an  analy¬ 
tical  solution  for  film  condensation  on  a  horizontal  plate  at 
a  constant  surface  temperature.  Their  analysis  and  resulting 
differential  equation  is  identical  to  the  development  in  Chapter 
II  for  a  smooth  cylindrical  condenser  if  the  acceleration  due 
to  gravity  is  replaced  by  a  radial  acceleration  term.  This 
being  the  case,  this  modification  was  made  and  the  analytical 
solution  was  used  as  a  reference  for  comparison. 

For  the  test  runs  of  the  finite  element  solution,  a  constant 
surface  temperature  was  assumed.  In  addition,  a  value  for  the 

maximum  film  thickness  (6*  „)  and  minimum  film  thickness  (&  •  ) 

max  min 

at  the  overfall  were  required.  The  value  for  6*  was  determined 

max 

based  on  a  relationship  developed  by  Leppert  and  Nimmo  [Ref.  8]. 
The  value  for  6m^n  was  arbitrarily  chosen. 

For  identical  geometry,  surface  temperature  and  maximum  and 
minimum  film  thicknesses,  the  results  of  both  analyses  were 
identical.  In  order  to  develop  confidence  in  the  finite  element 

A  ^  A 

solution,  <5_.  was  varied  from  0.10  •  to  0.97  •  <5 _ .  The 

min  max  max 

resulting  profiles  agreed  at  all  locations  along  the  length  of 
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the  condenser.  However,  when  the  temperature  was  permitted  to 
vary  along  the  length  of  the  condenser  in  the  finite  element 
solution  and  a  corresponding  average  temperature  was  used  in  the 
analytical  solution  of  Leppert  and  Nimmo  [Refs.  8  and  9],  the 
profiles  were  no  longer  in  agreement.  This  was  to  be  expected, 
particularly  in  the  case  where  there  was  a  sixteen  degree 
Fahrenheit  variation  along  the  length  of  the  condenser.  This 
substantial  temperature  variation  resulted  in  a  significant 
variation  of  fluid  properties  which  would  account  for  a  differ¬ 
ence  in  film  profile. 

Due  to  the  agreement  between  the  finite  element  solution 
and  the  analytical  solution  of  Leppert  and  Nimmo  [Refs.  8  and 
9]  for  a  constant  surface  temperature,  it  was  decided  that  the 
finite  element  solution  does  provide  a  satisfactory  representa¬ 
tion  of  the  film  profile.  This  being  the  case,  the  finite 
element  solution  was  then  incorporated  into  the  code  to  provide 
the  film  profile  for  the  cylindrical  condenser. 

Once  the  finite  element  solution  of  the  film  profile  was 
verified,  the  heat  transfer  analysis  could  be  accomplished. 

The  analysis  considered  both  copper  and  stainless  steel  con¬ 
densers  with  the  following  four  geometries:  a)  tapered-smooth, 
b)  tapered-axially  finned,  c)  cylindrical-smooth  and  d)  cyl- 
indrical-axially  finned.  Table  I  lists  the  geometric  para¬ 
meters  held  constant  for  all  analyses.  In  all  cases,  the 
working  fluid  was  water. 
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TABLE  I 


Condenser  Geometric  Parameters 
Held  Constant  During  All  Analyses 

condenser  length  =  8.500  inches 

minimum  radius  =  0.51575  inches 

wall  thickness  =  0.05118  inches 

In  addition,  the  following  geometric  parameters  were  also 
utilized  when  required.  This  requirement  was  based  on  the 
condenser  geometry  being  considered,  i.e.,  tapered-axially 
finned. 


TABLE  II 

Condenser  Geometric  Parameters 
Applied  as  Required 

height  of  fin  -  0.05906  inches 

fin  half  angle  =  26.565  degrees 

condenser  cone 

half  angle  =  1.00  degrees 

In  the  analysis,  the  heat  transfer  rate  was  determined  for 
the  four  different  geometries  listed  above  for  both  copper  and 
stainless  steel.  The  ambient  temperature  was  set  at  60.0°F 
and  the  heat  transfer  rate  was  determined  for  each  possible 
combination  of  the  operating  parameters  given  in  Table  III. 


62 


TABLE  III 


Operating  Parameter  Matrix 


Rotational 

Speed 

(RPM) 

Heat  Transfer 
Coefficient 
(btu/hr-ft2 -F) 

Saturation 
Temperature 
(degree  F) 

700.0 

100.0 

90.0 

1400.0 

500.0 

120.0 

2800.0 

1000.0 

150.0 

180.0 

Thus,  for  each  condenser  geometry,  there  was  a  total  of  72 
analyses,  36  for  the  case  of  the  condenser  with  a  copper  wall 
and  36  for  the  case  of  the  condenser  with  a  stainless  steel 
wall . 

The  first  condenser  geometry  considered  was  a  smooth  con¬ 
denser.  Figures  8-13  compare  the  heat  transfer  rates  of  smooth 
cylindrical  condensers  with  those  of  smooth  tapered  condenser. 

In  particular,  Figures  8,9,  and  10  indicate  the  results  of  the 
analyses  of  smooth  copper  condensers  at  rotational  speeds  of 
700,  1400  and  2800  revolutions  per  minute (RPM)  respectively. 
Figures  11,  12,  and  13  are  for  smooth  stainless  steel  conden¬ 
sers  at  700,  1400,  and  2800  RPM  respectively.  For  both  stain¬ 
less  steel  and  copper  smooth  condensers,  the  following  general 
observations  apply:  a)  For  the  same  external  heat  transfer 
coefficient,  the  heat  transfer  rate  for  the  cylindrical  smooth 
condenser  is  less  than  the  equivalent  tapered  condenser,  b)  As 
the  external  heat  transfer  coefficient  increases,  this  difference 


in  heat  transfer  rate  also  increases,  c)  The  rotational  speed 
has  a  greater  effect  on  the  tapered  condenser  heat  transfer 
rate.  For  example,  the  maximum  heat  transfer  rate  of  a  copper 
tapered  condenser  will  increase  by  a  factor  of  1.67  when  the 
rotational  speed  is  increased  from  700  to  2800  RPM.  In  the 
cylindrical  copper  condenser,  for  the  same  change  in  rotational 
speed,  the  heat  transfer  rate  only  increases  by  a  factor  of  1.51. 

These  same  observations  hold  true  for  the  smooth  stainless 
steel  condensers.  But,  due  to  a  greater  thermal  resistance  in 
the  wall  of  the  stainless  steel  condenser,  the  heat  transfer 
rates  are  less  for  all  cases  considered.  It  should  be  noted 
that  the  thermal  conductivity  of  stainless  steel  is  only  4%  the 
thermal  conductivity  of  copper.  This  accounts  for  the  increased 
thermal  resistance. 

Figures  14,  15,  and  16  compare  axially  finned  cylindrical 
with  axially  finned  tapered  copper  condensers  at  700,  1400  and 
2800  RPM  respectively.  Note  that  the  heat  transfer  rates  of 
the  cylindrical  condensers  are  only  slightly  less  than  those  of 
the  tapered  condensers.  This  is  because  the  heat  is  primarily 
transferred  through  the  extended  surface,  i.e.,  the  fin.  Thus 
the  film  condensate  in  the  trough  has  less  effect  on  the  heat 
transfer  rate.  In  fact,  the  average  difference  in  heat  transfer 
rate  for  700  RPM  and  an  external  heat  transfer  coefficient  of 
100  Btu/hr-ft2-F  is  12.65%.  This  difference  increases  to  15.6% 
as  the  heat  transfer  coefficient  is  increased  to  1000  Btu/hr- 
ft2-F.  However,  as  the  rotational  speed  is  increased  to  2800 
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RPM,  the  corresponding  average  differences  in  heat  transfer 
rates  decrease  to  12 . S 8 %  and  13.331  respectively. 

Figures  17,  18,  and  19  compare  the  heat  transfer  rates  of 
axially  finned  stailness  steel  cylindrical  condensers  with 
those  of  axially  finned  stainless  steel  tapered  condensers 
at  700,  1400  and  2800  RPM  respectively.  Note  the  difference 
in  heat  transfer  rate  increases  as  the  external  heat  transfer 
coefficient  increases  more  than  in  the  case  of  the  copper 
condensers.  At  the  low  heat  transfer  coefficient,  the  limit¬ 
ing  thermal  resistance  is  that  of  the  external  surface  (1/h  t)  . 
However,  as  the  external  heat  transfer  coefficient  is  increased, 
the  limiting  thermal  resistance  becomes  that  of  the  wall  due 
to  the  low  thermal  conductivity  of  stainless  steel. 

Another  observation  to  be  noted  is  the  fact  that  the  rota¬ 
tional  speed  has  a  greater  effect  on  the  heat  transfer  rate  of 
the  cylindrical  condensers  than  on  the  tapered  condensers.  As 
the  rotational  speed  increases,  the  film  thickness  decreases 
due  to  the  greater  centrifugal  force  exerted  on  the  film.  For 
the  cylindrical  condenser,  the  maximum  film  thickness  is  much 
greater  than  for  a  tapered  condenser.  For  example,  an  axially 
finned  stainless  steel  cylindrical  condenser  rotating  at  1400 
RPM  with  an  external  heat  transfer  coefficient  of  1000  Btu/hr- 
ft2-F  has  a  maximum  film  thickness  twice  that  of  a  corresponding 
tapered  condenser.  This  being  the  case,  higher  rotational 
speeds  will  have  a  greater  effect  on  the  greater  film  thickness 
and  the  difference  in  heat  transfer  rates  will  decrease. 
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Figures  20,  21,  and  22  compare  the  heat  transfer  rates  of 
copper  smooth  cylindrical  condensers  with  copper  axially  finned 
cylindrical  condensers  at  700,  1400  and  2800  RPM  respectively. 
The  figures  indicate  that  for  low  external  heat  transfer  coef¬ 
ficients,  little  is  to  be  gained  by  the  addition  of  axial  fins. 
But,  as  the  external  heat  transfer  coefficient  increases,  the 
advantage  becomes  significant.  As  an  example,  consider  Figure 
21.  For  h=100  Btu/hr- ft2 - F ,  axial  finning  increase  the  heat 
transfer  rate  by  161.  But,  for  an  external  heat  transfer  coef¬ 
ficient  of  1000  Btu/hr- ft 2 - F ,  the  heat  transfer  rate  increases 
by  194%. 

Note  also,  that  as  rotational  speed  increases  the  advan¬ 
tage  to  be  gained  by  axial  finning  decreases.  This  can  be 
explained  by  the  fact  that  for  an  axially  finned  condenser, 
the  majority  of  the  heat  is  transferred  by  the  fin  surface. 

As  the  rotational  speed  increases,  the  film  thickness  in  the 
trough  decreases.  This  in  turn  will  expose  slightly  more  fin 
surface  area.  On  the  other  hand,  for  the  smooth  condenser, 
the  decrease  in  the  film  thickness  will  have  a  greater  effect 
on  heat  transfer  rate  in  that  the  thermal  resistance  of  the 
film  has  decreased.  Comparing  the  two  geometries,  the  change 
in  overall  thermal  resistances  for  the  smooth  condenser  will 
be  greater  than  that  for  the  axially  finned  condenser  account¬ 
ing  for  the  slight  decrease  in  advantage  with  increasing 
rotational  speed. 
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Figures  23,  24,  and  25  correspond  to  Figures  20,  21,  and 
22  but  for  stainless  steel  condensers.  The  results  are  similar 
to  the  copper  situation  discussed  above,  by  the  effect  caused 
by  the  increasing  external  heat  transfer  coefficient  is  not 
as  dominant  due. to  the  high  thermal  resistance  of  the  stainless 
steel  wall  material. 

Figure  26  compares  a  smooth  cylindrical  copper  condenser 
with  a  smooth  cylindrical  stainless  steel  condenser  at  1400 
RPM.  As  to  be  expected,  the  heat  transfer  rate  of  the  copper 
condenser  is  greater  than  that  of  the  stainless  steel  condenser 
due  to  the  difference  in  the  thermal  conductivity  of  the  two 
materials.  The  difference  in  the  heat  transfer  rate  is  least 
for  a  low  external  heat  transfer  coefficient  where  the  external 
thermal  resistance  is  dominant.  As  the  heat  transfer  coeffi¬ 
cient  increases,  the  thermal  resistance  of  the  wall  of  the 
condenser  becomes  more  important  resulting  in  an  increasing 
difference  between  the  two  condensers. 

Figure  27  provides  a  comparison  of  axially  finned  cylindrical 
copper  and  stainless  steel  condensers  at  1400  RPM.  The  advantage 
of  copper  over  stainless  steel  is  obvious.  Note  that  the  heat 
transfer  rate  for  the  copper  condenser  at  500  Btu/hr-ft2-F  is 
nearly  identical  to  the  stainless  steel  condenser  heat  transfer 
rate  at  1000  Btu/hr/-ft2 -F  indicating  the  advantage  of  copper 
over  stainless  steel. 

The  final  analysis  compared  the  heat  transfer  rate  of  an 
axially  finned  copper  condenser  with  a  triangular  fin  profile 


to  the  heat  transfer  rate  of  an  axially  finned  condenser  with 
a  rectangular  fin  profile.  The  rectangular  fin  was  assumed  to 
have  an  adiabatic  tip.  This  comparison  was  accomplished  for 
both  cylindrical  and  tapered  condensers.  Table  IV  lists  the 
parameters  used  in  the  analysis.  These  parameters  are  in  addi 
tion  to  those  listed  in  Table  I.  This  comparison  was  only 
accomplished  for  the  one  set  of  operating  parameters  listed  in 
Table  IV.  Note  that  the  operating  parameters  chosen  were  the 
median  values. 


TABLE  IV 

List  of  Parameters  Used  in 
Rectangular/Triangular  Fin  Profile  Comparison 


Heat  Transfer  Coefficient 
Rotational  Speed  = 

Saturation  Temperature  = 

Fin  Height  = 

Table  V  lists  the  results  of  the  analyses 


500  Btu/hr-ft2-F 
1400  RPM 
120.0  degrees  F 
0.05906  inches 


TABLE  V 


Results  of  Rectangular/Triangular 
Fin  Profile  Comparison 


CONDENSER  GEOMETRY  FIN  PROFILE 


tapered  triangle 

rectangle 
triangle 


Q  (Btu/hr) 
6168.93 
6143.03 
5358.4 


cylindrical 


Note  that  the  heat  transfer  rate  varied  by  only  0.4%  for  both 
the  tapered  and  cylindrical  condensers.  The  reason  for  this 
insignificant  variation  is  provided  in  Table  VI  and  VII. 

Table  VI  lists  the  convective  surface  temperature  distribution 
of  a  symmetric  unit  section  of  a  tapered  condenser  at  the  middle 
increment  of  the  condenser.  Temperature  location  1  is  located 
at  the  tip  of  the  fin.  For  the  rectangular  profile,  temperature 
location  1  is  located  at  the  intersection  of  the  adiabatic 
surface  (the  tip  of  the  fin) ,  and  the  vertical  surface  of  the 
fin.  Temperature  location  5  is  located  at  the  base  of  the  fin. 
Locations  6-9  are  located  in  the  trough  region  and  the  remaining 
temperature  locations  are  along  the  external  surface  of  the 
section.  Thus  temperatures  1-5  provide  the  temperature  distri¬ 
bution  along  the  convective  surface  of  the  fin.  Note  that  the 
temperature  distribution  is  lower  for  the  rectangular  fin. 

Thus,  the  driving  force  for  heat  transfer,  the  temperature 
difference  between  the  saturation  temperature  and  the  surface 
temperature  is  greater  for  the  rectangular  fin.  Thus,  in  spite 
of  the  fact  that  the  fin  surface  area  for  heat  transfer  has  de¬ 
creased  by  11%,  the  average  surface  temperature  difference  has 
increased  by  25%.  It  should  also  be  noted  that  the  average 
heat  transfer  coefficient  for  the  fin  surface  for  the  rectan¬ 
gular  fin  also  decreases  by  12%.  However,  the  increase  in  temper 


ature  difference  is  the  dominant  change  that  only  allows  a  0.4% 
decrease  in  heat  transfer  rate. 


TABLE  VI 


Surface  Temperature  Distribution  for  Rectangular 
and  Triangular  Axially  Finned  Tapered  Copper 
Condensers  at  Middle  Increment  of  Condensers 


Temperature 

Location 


Triangular  Fin 
Temperature  (F) 

118.85 

118.45 

118.04 


Rectangular  Fin 
Temperature  (F) 

117.95 

117.88 

117.68 


117.60 
117.05 
116 . 82 

116.61 
116.59 
116.48 
116.47 

116.46 

116.46 
116.44 
116.42 
116.31 
116.31 
116.27 


116.26 


117.36 


116.59 

116.45 

116.35 

116.35 

116.54 

116.53 

116.23 

116.21 

116.20 

116.17 

116.12 

116.07 

116.04 

116.03 


TABLE  VII 

Surface  Temperature  Distribution  for  Rectangular 
and  Triangular  Axially  Finned  Cylindrical  Copper 
Condensers  at  Middle  Increment  of  Condensers 


Temperature 

Location 


Triangular  Fin 
Temperature  (F) 

118.55 
118.09 
117.62 
117.10 

116.55 
116.38 
116.28 
116.22 
116.20 
116.03 
116.02 
116.01 
115.99 
115.97 
115.94 
115.91 
115.88 
115.88 


Rectangular  Fin 
Temperature  (F) 

117.55 

117.44 

117.25 

116.87 

116.36 

116.18 

116.07 

116.01 

115.99 

115.80 

115.79 


115.79 


115.77 


115.75 

115.72 

115.67 

115.67 


115.66 


Table  VII  indicates  a  similar  situation  exists  in  the 
cylindrical  condenser.  Again  note  the  decrease  in  the  average 
surface  temperature  of  the  fin  for  the  rectangular  profile. 
The  surface  area  and  heat  transfer  coefficient  has  decreased, 
but  the  increased  temperature  difference  compensates  for  these 
changes,  limiting  the  overall  decrease  in  heat  transfer  rate. 


Heat  Transfer  Rate  (Btu/hr) 


Saturation  Temperature  (°F) 


1. 

^ext 

2. 

^ext 

3. 

h 

ext 

100  Btu/hr-ft2 -F 
500  Btu/hr-ft2 -F 
1000  Btu/hr-ft2 -F 


Figure  8.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Smooth  Cylindrical  and  Tapered  Copper  Con¬ 
densers  at  700  RPM 
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Figure  9.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Smooth  Cylindrical  and  .Tapered  Copper  Con¬ 
densers  at  1400  RPM 
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Figure  11.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Smooth  Cylindrical  and  Tapered  Stainless 
Steel  Condensers  at  700  RPM 
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Figure  12.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Smooth  Cylindrical  and  Tapered  Stainless 
Steel  Condensers  at  1400  RPM 
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Figure  13.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Smooth  Tapered  and  Cylindrical  Stainless 
Steel  Condensers  at  2800  RPM 
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Figure  14.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Axially  Finned  Cylindrical  and  Tapered  Copper 
Condensers  at  700  RPM 
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Figure  15.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Axially  Finned  Cylindrical  and  Tapered  Copper 
Condensers  at  1400  RPM 
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Figure  16.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Axially  Finned  Cylindrical  and  Tapered  Copper 
Condensers  at  2800  RPM 
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7.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Axially  Finned  Cylindrical  and  Tapered  Stain¬ 
less  Steel  Condensers  at  700  RPM 
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Figure  18.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Axially  Finned  Cylindrical  and  Tapered  Stain¬ 
less  Steel  Condensers  at  1400  RPM 
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Figure  19.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Axially  Finned  Cylindrical  and  Tapered  Stain 
less  Steel  Condensers  at  2800  RPM 
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Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Smooth  and  Axially  rinned  Copper  Cylindrical 
Condensers  at  700  RPM 
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Figure  21.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Smooth  and  Axially  Finned  Copper  Cylindrical 
Condensers  at  1400  RPM 
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Figure  22.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Smooth  and  Axially  Finned  Copper  Cylindrical 
Condensers  at  2800  RPM 
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Figure  23.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Smooth  and  Axially  Finned  Stainless  Steel 
Cylindrical  Condensers  at  700  RPM 
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Figure  24.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Smooth  and  Axially  Finned  Stainless  Steel 
Cylindrical  Condensers  at  1400  RPM 
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Figure  25.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Smooth  and  Axially  Finned  Stainless  Steel 
Cylindrical  Condensers  at  2800  RPM 
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Figure  26.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Copper  and  Stainless  Steel  Smooth  Cylindrical 
Condensers  at  1400  RPM 
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7.  Heat  Transfer  Rate  versus  Saturation  Temperature 
for  Copper  and  Stainless  Steel  Axially  Finned 
Cylindrical  Condensers  at  1400  RPM 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  heat  transfer  rate  of  a  cylindrical  condenser  is  less 
than  an  equivalent  tapered  condenser.  This  decrease  in  heat 
transfer  rate  is  most  significant  in  a  smooth  condenser,  where, 
depending  on  the  external  heat  transfer  coefficient  and  rota¬ 
tional  speed,  can  be  as  great  as  451.  This  decrease  in  heat 
transfer  rate  becomes  less  significant  for  an  axially  finned 
condenser  where  the  average  decrease,  for  the  range  of  heat 
transfer  coefficients  examined,  was  13%  for  a  copper  condenser 
and  18%  for  a  stainless  steel  condenser.  When  such  factors 
as  the  cost  and  difficulty  in  manufacturing  tapered  axially 
finned  condensers  are  considered,  this  13%  decrease  in  heat 
transfer  rate  becomes  tolerable.  From  a  practical  standpoint, 
development  and  analysis  of  cylindrical  axially  finned  con¬ 
densers  should  be  encouraged  by  these  results.  If  environ¬ 
mental  conditions  permit,  copper  should  be  preferred  over 
stainless  steel  due  to  its  exceptionally  high  thermal  con¬ 
ductivity  and  resulting  higher  heat  transfer  rate. 

B.  RECOMMENDATIONS 

1)  Build  and  experimentally  test  both  smooth  and. axially 
finned  cylindrical  condensers  to  obtain  experimental  data 
for  comparison  with  results  of  this  analysis. 
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2)  Develop  models  for  rectangular  and  trapezoidal  fin 
profiles  with  non- adiabatic  tips  and  incorporate  into  code 

3)  Experimentally  test  axially  finned  cylindrical  con 
densers  with  rectangular  and  trapezoidal  fin  profiles  to 
obtain  data  for  comparison  with  theoretical  results. 


APPENDIX  A 


FILM  PROriLE  FINITE  ELEMENT  SOLUTION 

A.  SMOOTH  CONDENSER 

The  analysis  of  the  film  profile  in  a  smooth  cylindrical 
condenser  developed  in  Chapter  II  resulted  in  the  following 
ordinary,  nonlinear,  second  order,  differential  equation: 


r  *  _fL  *  *  3  1 

5  37  L33T  5  J 


5kf(Tsat;  -  Tw)  u 

Pr2’w2r 

f  fg 


(eqn  A.l) 


This  equation  can  be  rearranged  and  expanded  to  yield: 


*4  d2,s* 


il  +  35  * 3  (^0  =  -K 

dx2  dx 


(eqn  A. 2) 


where  K  *  3kf'^TSat  "  Tu^w 


p <:2(u2 r  E 


fg 


The  statement  of  the  problem  for  the  formulation  of  the 
Finite  Element  Method  is : 


5*“  +  35  * 3  (f^)  =  -K  (eqn  A.  3) 

dx2  dx 

with  the  following  boundary  conditions: 

a)  at  x  =  0,  i*  =  x* 

b)  at  x  =  0,  d§*/dx  =  0 
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i  RESOLUTION  TEST  CHART 

BUREAU  OF  STANDARDS-  1963-A 


The  domain  (length  of  the  condenser)  is  divided  into  ele¬ 
ments  of  length  Ax  with  the  exception  of  the  first  and  final 
elements  which  have  length  Ax/2,  where  Ax  is  the  length  of  the 
condenser  divided  by  the  number  of  axial  increments  (NDIV) . 

A  system  nodal  point  is  located  at  each  end  of  an  element. 
Thus,  system  nodal  point  1  is  located  at  x=0  and  the  last  sys¬ 
tem  nodal  point,  which  is  equal  to  the  number  of  elements  plus 
1  is  located  at  the  overfall  into  the  evaporator,  x=L.  All 
internal  nodal  points  are  located  at  a  position  corresponding 
to  the  midpoint  of  the  axial  increment. 

Define  the  approximate  value  of  the  film  thickness  in  the 
following  manner: 


7*6 


* 

n 


(x) 


h  m 

r  g. d.  =  g  d 

1  1  X 


(eqn  A.  4) 


where  7  =  the  approximate  value  of  the  film  thickness  (6*) . 
=  the  global  basis  functions, 
n  =  the  number  of  system  nodal  points, 
d  =  the  solution  vector. 

On  an  element  level,  equation  (A. 4)  becomes: 


4 

"5-  ~  ofx)  =  Z  g.d,  (eqn  A. 5) 

e  e  i=l  1  xe 


where  g^  =  the  local  basis  functions. 
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Define  two  degrees  of  freedom  at  each  nodal  point,  i.e., 
both  "S’  and  dS/dx  are  continuous.  Thus,  the  local  basis  func¬ 
tions  used  in  the  finite  element  solution  are: 


111  +  ill 

Z2  z3 


(eqn  A. 6) 


(eqn  A. 7) 


(eqn  A.  8) 


(eqn  A. 9) 


where  Z  =  the  length  of  an  element. 

gj^  and  g3  =  magnitude  basis  functions. 
g2  and  g4  =  slope  basis  functions. 

Substituting  the  approximate  film  thickness  into  the 
differential  equation  (A. 3)  results  in: 


j'  =  -K„  (eqn  A. 10) 

e  dx  e  dx  e 

To  remove  the  "nonlinearity"  from  the  problem,  equation  (A. 10) 
is  modified  in  the  following  manner: 


d2T  dT 

n4 - —  +  3n3r)f  -  =  -K  (eqn  A.  11) 

dx2  dx  e 


9-7 


n  is  defined  as  the  approximate  value  of  the  film  thickness 
from  the  previous  iteration.  In  like  manner,  n '  is  defined 
as  the  approximate  value  of  the  rate  of  change  of  film  thick¬ 
ness  with  respect  to  x  from  the  previous  iteration. 

Forming  the  residual  of  equation  (A. 11)  yields: 

Rn  =  n**  (GT"d)  +  3nV(GT  d)  +  K  (eqn  A.  12) 

Invoking  the  Galerkin  criterion  for  the  determination  of  the 
solution  vector  d,  i.e. 

/G^Rndx  =  C  i  =  l,2,3,...n  (eqn  A. 13) 

yields : 

nVGGT  ddx  +  3nV/  GG  ddx  +  K/Gdx  =  0  (eqn  A. 14) 


Each  of  the  integrals  in  equation  (A. 14)  are  defined  in  the 
following  manner: 


T"  n  T" 

/GG  ddx  =  l  gg  d  dx 
1  e 


(eqn  A. 15) 


Or  on  an  elemental  level: 


??  _  ff  _  tT  _  If 


/ggT  dedx  =  / 


r«ll  <gl  g2  §3  §3  " 
g2 
g3 
g4 


dgdx  (eqn  A. 16) 


This  integration  results  in  the  following  4x4  local  elemental 
A  matrix  for  any  element: 
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manner , 

let 

/GGr  ddx  =  Z  fzggT'  dedx 


(eqn  A. 17) 


and  on  an  elemental  level 


fzgg  d  dx  =  fi 


fgl 

g2 

g3 

g4 


<gl  g2  g3  g4  > 


d  dx 
e 


Ceqn  A. 18) 


This  integration  results  in  the  following  4x4  local  elemental 
B  matrix  for  any  element: 
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Lastly,  let 


n 


/  Gdx  =  2  f  gdx 

1  *  ~ 


(eqn  -A.  19) 


or,  on  an  elemental  level  this  becomes 


/£gdx  -  /j 


*1 

g2 

g3 

g4 


dx 


(eqn  A. 20) 


This  integration  results  in  the  following  local  elemental 
column  F  vector: 


l 

1 


[FI. 


12 


h 


-r 

n 


Thus,  for  a  given  element,  equation  (A. 12)  becomes 


[V[A]e  +  3  3  ne'  [B]e]  ^e  =  (eqn  A*21) 

As  mentioned  above,  n  and  n'  ate  the  approximate  values  of  5* 

and  dfi*/dx  respectively  from  the  previous  iteration.  For  the 

initial  iteration,  n  is  set  equal  to  5*  and  n '  is  set  equal 

max 

to  0 . 

Each  elemental  matrix  is  placed  in  a  global  system  [A] 


matrix  with  the  location  in  the  matrix  based  upon  the  local/global 


nodal  point  correspondence.  For  example,  for  element  Nr.  2 
the  nodal  points  corresponding  to  local  nodal  points  1,2,3, 
and  4  are  global  nodal  points  3,4,5,  and  6  respectively. 

Therefore,  the  sum  of  the  two  4x4  elemental  matrices  ([A]  and 
[B])  and  multiplying  constants  will  form  a  single  4x4  local 
matrix.  Element  (1,1)  of  this  local  matrix  will  be  placed  in 
the  global  A  matrix  location  (3,3)  and  element  (4,4)  of  the  local 
elemental  matrix  will  be  placed  in  the  global  A  matrix  location 
(6,6).  Ultimately,  the  following  global  system  would  be  assembled 

[A)mXm  •  d  -  l  (eqn  A. 22) 

Note  that  the  global  A  matrix  would  be  an  mxm  size  matrix 
where  m  is  equal  to  twice  the  number  of  system  nodal  points  to 
account  for  the  two  degrees  of  freedom  at  each  nodal  point. 

In  a  similar  fashion,  d  and  F  would  be  column  vectors  of  size  m. 

Once  the  system  is  assembled  the  boundary  conditions  are 
applied.  This  is  done  in  the  following  manner:  For  boundary 
condition  (a),  A(l,l)  is  set  equal  to  1.0  and  the  remaining 
elements  in  the  first  row,  i.e.  A  (l,i)  i=2,3,4,...m  are  set 
equal  to  0.0.  Then  F(l)  is  set  equal  to  6*  •  6 * _ „  is  ini- 

tially  determined  by  a  relationship  developed  by  Leppert  and 
Nimmo  [Refs'.  8  and  9],  This  establishes  the  value  of  d(l)  as 
<5*  .  In  a  similar  manner,  boundary  condition  (b)  is  applied 

by  setting  A(2,2)  equal  to  1.0  and  all  other  second  row  elements 
of  the  global  A  matrix  are  set  equal  to  0.0.  Then,  F(2)  is  set 
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equal  to  0.0.  This  establishes  the  solution  vector  d(2)  as  0.0. 
Note  that  d(2)  corresponds  to  the  slope  at  the  first  nodal  point. 

One  additional  boundary  condition  is  required;  this  is  the 
value  of  the  film  thickness  at  x=L.  This  boundary  condition  is 
necessary  to  completely  specify  the  problem.  The  value  of  the 
film  thickness  at  the  overfall  may  take  on  any  value  depending 
on  the  geometry  at  the  overfall.  In  the  case  of  this  analysis, 
this  value  was  taken  as  0.25  •  6*ax-  This  value  was  chosen  for 
the  following  reason:  Leppert  and  Nimmo  [Refs.  8  and  9],  in 
a  similar  analysis  for  laminar  film  condensation  on  a  horizontal 
surface  derived  an  analytical  solution  to  equation  (A.l),  assum¬ 
ing  a  constant  surface  temperature.  They  found  the  film  profiles 
for  6*  with  the  overfall  value  less  than  0.40  .  5*ax  were  essen¬ 
tially  constant  and  thus  any  value  of  6*  at  the  overfall  less 
than  0.40  •  6*  would  result  in  the  same  profile.  In  the  veri- 
fication  of  the  finite  element  solution,  not  only  was  this  found 
to  be  the  case,  but  it  was  also  found  that  the  heat  transfer 
rate  was  relatively  insensitive  to  the  shape  of  the  film  profile. 
In  fact,  it  was  found  that  the  film  thickness  at  the  overfall 
could  be  increased  to  a  value  as  great  a  0.90  .  6*ax  and  the 
resulting  variation  in  heat  transfer  rate  was  only  10%.  This 


being  the  case,  the  value  of  0.25  •  <S  *ax  was  arbitrarily  chosen 

for  <$  •  _  . 
mxn 

The  third  boundary  condition  was  applied  by  setting  A(m-l,i)= 
0.0  where  i=l,2,3...m.  Then,  the  global  matrix  element  A(m-l,m-l) 
was  set  equal  to  1.0.  Finally,  F(m-l)  was  set  equal  to  0.25*6*  • 
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Once  all  three  boundary  conditions  were  applied,  the  sys¬ 
tem  given  by  equation  (A. 22)  was  solved  for  d.  The  values  of 
the  approximate  film  thickness,  i.e.,  d(i)  ,  i=l , 3 , 5 , . . .m- 1  are 
then  compared  to  the  values  of  film  thickness  from  the  previous 
iteration.  If  the  relative  difference  is  less  than  or  equal 
to  a  specified  convergence  criterion  (i.e.,  0.0001)  at  all 
nodal  points,  convergence  is  met  and  the  latest  values  of  d 
are  the  solution  values  of  5*. 

If  convergence  is  not  met,  the  values  of  d  are  saved  for 
the  next  iteration  where  they  are  used  to  determine  n  and  n' 
as  discussed  above.  This  iterative  process  is  continued  until 
convergence  is  met  or  until  a  maximum  number  of  iterations  have 
occurred. 

B.  AXIALLY  FINNED  CONDENSER 

The  finite  element  solution  for  the  film  profile  in  an 
axially  finned  cylindrical  condenser  is  very  similar  to  that 
of  a  smooth  cylindrical  condenser.  For  this  reason,  only  the 
variations  in  the  development  will  be  addressed.  From  Chapter 
II,  the  differential  equation  for  mass  flow  rate  in  an  axially 
finned  cylindrical  condenser  is  given  by: 


<5*  [^-  (eS*3+  5*1* tana)  ] 


3kf (Tsat  ~  Tw)  yte 
°f2(i)2r  ^fg 


-25*cosa 


4kfCTsat  '  Tav,>z 


l  "f 


pjr2o)2r  h 


fg 


cosa 


(eqn  A. 23) 
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This  equation  can  be  rearranged  and  expanded  to  yield: 


^  (€«*•♦  S*stana)  +  (3eJ*3  ^  *  46*'taim 


-  -Kj  -K2«* 


(eqn  A. 24) 


where  K,  =  3kftTsat  '  Tw)u£ 


Pf2»!r  K 


*  2cosa 


4kf(Tsat  -  Tavg)"zl 
p^2ui2r  cosot 


Substituting  equation  (A. 5)  into  equation  (A. 24)  results  in 


d27  _4  _s  d6  _3  d6  4  d3\ 

(eSe  *  4e  tana)  *  ~  (3e“e  37"  *  4Te  tan“  3T"5 


K1  -  K2Se 


(eqn  A. 25) 


To  remove  the  "nonlinearity"  from  the  problem,  equation  (A. 25) 
is  modified  in  the  following  manner" 


d2?  dff 

-  (eY1*  +Y 5 tana)  +  -  (3cy3Y,+  4Y4Y'tana) 


K,  -  K^Y 


(eqn  A. 26) 


it 


Here,  y  and  y'  are  defined  by  the  following  relationships: 


Y  =  +  R*C^  (eqn  A-27) 

y *  *  "5\  +  R*f5^  _^)  (eqn  A. 28) 


where  =  approximate  value  of  the  film  thickness  for 

the  present  iteration. 

approximate  value  of  film  thickness  from  the 
previous  iteration. 

_ 1 

6^  =  approximate  value  of  the  derivative  of  the  film 
thickness  for  the  present  iteration 
approximate  value  of  the  derivative  of  the  film 
thickness  from  the  previous  iteration. 

R  *  relaxation  factor. 

These  two  variables,  y  and  y',  are  in  actuality,  adjusted  approx 
imation  of  film  thickness  and  derivative  of  film  thickness  re¬ 
spectively.  This  adjustment  is  required  in  order  to  converge  to 
a  solution. 

The  finite  element  solution  is  now  identical  to  that  of  the 
smooth  condenser,  that  is,  the  residual  is  formed,  the  Galerkin 
criterion  is  invoked,  and  identical  4x4  local  matrices  are 
derived.  Finally,  the  equivalent  of  equation  (A. 21)  is  formed. 


[(eye4  +  Ye5tana)[A]e  +  (3eye  y^  +  Ye4Yg  tana)  [B]g]  •  dg 


*  <  -  Kle  *  K2e  V  fe 


(eqn  A. 29) 
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Notice  that  this  equation  has  two  forcing  terms.  The  additional 
term  (^Y)  resulted  from  the  nonlinearization  of  the  problem. 

Just  an  in  the  smooth  condenser  film  profile  solution,  the 
global  system  given  by  equation  (A. 22)  is  formed,  the  boundary 
conditions  applied  and  the  system  solved  for  a  solution  vector 
d.  The  iterative  process  is  continued  until  convergence  is  met. 
With  each  iteration,  y  and  y'  are  updated  and  used  for  the 
next  iteration.  When  convergence  is  met,  the  latest  values  of 
dfi)  ,  i=l , 3 , 5 , . . .m- 1  are  the  solution  values  of  6*. 


APPENDIX  B 


USER'S  MANUAL 

This  appendix  describes  the  data  cards  required  to  use 
the  computer  code. 

The  data  is  divided  into  "blocks"  for  convenience.  Each 
page  of  this  user's  guide  is  a  separate  block.  For  each  block, 
a  general  description,  the  required  format,  and  appropriate 
comments  are  provided. 

It  is  imperative  that  input  data  be  consistent  or  errors 
will  result.  For  example,  if  a  smooth  geometry  is  being  analyzed, 
the  finite  element  parameters  must  also  result  in  a  smooth  model. 
In  addition,  all  data  fields  must  be  filled  with  an  input  value, 
even  if  that  value  is  not  needed  for  the  analysis.  For  example, 
in  a  smooth  analysis,  no  fin  half  angle  is  required  for  the  cal¬ 
culations;  however,  a  value  of  the  correct  format  must  be  pro¬ 
vided  in  the  fin  half  angle  field  or  an  INPUT/OUTPUT  error  will 
result . 
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DATA  BLOCK 


A 


DESCRIPTION:  FINITE  ELEMENT  PARAMETERS 
FORMAT :  8 1 S 

12345678 


NDIV  NCMREC  NCMTRI  NRIVFIN  NRWTRF  NCMTRF  NCOL  NPRNT 


FIELD  CONTENTS 

1  NDIV - Number  of  axial  increments.  Must  be  less  than 

or  equal  to  50. 

2  NCMREC-Number  of  columns  of  finite  elements  in  the  rec¬ 

tangular  section  of  fin.  May  be  equal  to  0  if 
fin  is  triangular  only. 

3  NCMTRI-Number  of  columns  of  finite  elements  in  tri¬ 

angular  section  of  fin.  May  be  equal  to  0  if 
only  rectangular  fin. 

4  NRWFIN-Number  of  rows  of  finite  elements  in  the  fin 

section  of  model.  May  be  equal  to  0  if  a  smooth 
condenser . 

Note:  If  triangular  or  trapezoidal  fin,  NRWFIN 

must  equal  NCMTRI. 

5  NRWTRF-Number  of  rows  of  finite  elements  in  wall  section 

of  the  model. 

6  NCMTRF-Number  of  columns  of  finite  elements  in  the 

trough  section  of  finned  model.  Set  equal  to 
0  if  a  smooth  condenser. 

7  NC0L---Total  number  of  columns  of  finite  elements. 

Must  be  equal  to  NCMREC+NCMTRI+NCMTRF  for  a 
finned  model. 

8  NPRNT- -Print  control  number. 

If  equal  to  1  correspondence  table  and  major 
elements  of  finite  element  model  will  be  printed 
out . 

If  equal  to  0,  output  will  be  suppressed. 
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DATA  BLOCK  B 


DESCRIPTION:  FLUID,  FIN  MATERIAL  SELECTION  PARAMETERS 

FORMAT :  215 

12345678 


DATA  BLOCK  C 


CLI  REAS El  THICKI  BFINI  CANGL  FNWTHI 


CONTENTS 

CLI - Condenser  length  (inches). 

RBASEI - Inside  radius  to  wall  of  condenser  at  condenser 
end  (inches). 

THICKI -Wall  thickness  (inches). 

BFINI --Fin  height  (inches).  Must  be  set  equal  to 
0.0  if  smooth  condenser. 


CANGL- - Condenser  half  angle  for  tapered  condenser 
(degrees).  Must  be  set  equal  to  0.0  for 
cylindrical  condenser. 

FNWTHI-Width  of  rectangular  portion  of  trapezoidal 
or  rectangular  fin  (inches).  If  triangular 
fin,  set  equal  to  0.0. 


DATA  BLOCK  D 


DESCRIPTION:  INTERNAL  FIN  GEOMETRY 

FORMAT:  2G10.5 

12345678 

FANGL  ETOEO 


FIELD  CONTENTS 

1  FANGL--Fin  half  angle  (degrees).  Set  equal  to  0.0 

for  smooth  condenser  or  rectangular  fin. 

2  ETOEO--Ratio  of  trough  width  to  fin  base  width. 

Determines  spacing  between  fins. 


Ill 


DATA  BLOCK  E 


DESCRIPTION:  CONVERGENCE  CRITERIA 
FORMAT:  2G10.5 

12345678 


CRIT  CRITDL 


FIELD  CONTENTS 

1  CRIT — Temperature  convergence  criterion.  Used  to 

determine  solution  of  two  dimensional  steady 
state  heat  conduction  problem. 

2  CRITDL-Mass  flow  convergence  criterion.  Used  only 

in  cylindrical  condenser  analysis  for  mass 
flow  convergence  test. 


r 
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DATA  BLOCK  F 


FIELD  CONTENTS 

1  RPM - Rotational  speed  (revolutions  per  minute). 

2  HINF - External  heat  transfer  coefficient 

(Btu/hr-ft  -deg  F)  . 

3  TINTL- -  Initial  temperature  estimate  (degrees  F) . 

4  TSAT - Saturation  temperature  (degrees  F) . 

5  TINF - Ambient  temperature  (degrees  F) . 


IUNITS  NFLAG1  NFLAG2  NFLAG3 


FIELD  CONTENTS 

1  IUNITS-Output  units  control  number. 

If  IUNITS  ■  0,  calculated  results  will  be  pro¬ 
vided  in  English  units. 

If  IUNITS  =  1,  input  parameters  will  be  repeated 
in  SI  units  and  calculated  results  will  be  pro¬ 
vided  in  SI  units. 

If  IUNITS  =  2,  input  parameters  will  be  repeated 
in  SI  units  and  output  results  will  be  provided 
in  both  English  and  SI  units. 

2  NFLAGl-The  first  axial  increment  at  which  the  parameters 

listed  under  remarks  will  be  provided  as  output. 

3  NFLAG2-The  final  increment  at  which  the  parameters 

listed  under  remarks  will  be  provided. 

4  NFLAG3-The  step  change  in  increments  between  NFLAG1. 

REMARKS 

1)  No  matter  what  value  of  IUNITS  is  used,  input  parameters 
will  always  first  appear  in  English  units. 

2)  For  increments  indicated  by  NFLAG  values,  the  following 
parameters  will  appear:  a)  x- coordinate ,  y-coordinate  and 
temperature  at  each  nodal  point,  b)  length  of  element,  heat 
transfer  coefficient  and  heat  rate  per  unit  length  for  each 
convective  boundary  element. 

3)  As  a  minimum,  the  values  of  1,  1,  1  must  be  provided  as 
input  for  NFLAG1,  NFLAG2 ,  and  NFLAG3  respectively. 
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NSOLVE 


FIELD  CONTENTS 

1  NSOLVE-Tapered  solution  control  number. 

For  tapered,  axially  finned  condenser,  set 
NSOLVE* 1 . 

For  tapered  smooth  condenser,  NSOLVE  must  be 
set  to  one  of  the  following  three  values: 

Set  NSOLVE  *  2  if  solution  of  film  thickness  is 
to  be  based  on  Ballback's  [Ref.  1] 
equation. 

Set  NSOLVE  =  3  if  solution  of  film  thickness  is 
to  be  based  on  Daniels  and  Al-Jumaily  [Ref.  13] 
equation,  neglecting  drag  terms. 

Set  NSOLVE  =  4  if  solution  of  film  thickness  is 
to  be  based  on  Daniel's  and  Al-Jumaily  [Ref. 

13]  equation,  with  drag  effects  included. 


REMARKS 


DATA  BLOCK  I 


DESCRIPTION:  CYLINDRICAL  CONDENSER  ANALYSIS  PARAMETERS 
FORMAT:  3I5,2GI5.10 

12345678 


NONCE  ITERMX  ITRPRT  RELAX  DELMAX 


FIELD  CONTENTS 

1  NONCE- -Single  iteration  parameter. 

If  NONCE  ■  1,  only  one  iteration  will  be 
permitted. 

If  NONCE  =  0,  iterations  will  continue  until 
convergence  or  maximum  number  of  iterations 
is  reached. 

2  ITERMX- Maximum  number  of  iterations  permitted  in  analysis. 

3  ITRPRT- Iteration  print  control  parameter. 

If  ITRPRT  =  1,  mass  flow  convergence  test  results 
will  be  provided  for  each  iteration. 

If  ITRPRT  =  0,  mass  flow  convergence  test  results 
will  only  be  provided  on  final  iteration. 

4  RELAX- -Relaxation  variable  used  in  finite  element  solution 

of  cylindrical  finned  film  profile. 

5  DELMAX- Initial  estimate  of  maximum  film  thickness  used  in 

solution  of  cylindrical  finned  film  profile. 

REMARKS 

1)  Above  parameters  are  only  used  in  cylindrical  condenser 
analysis . 

2)  Recommended  value  of  RELAX  is  0.80.  It  is  sometimes  necessary 
to  adjust  this  value  plus  or  minus  0.05  to  reach  film  profile 
convergence  at  small  film  thickness  values. 

3)  Input  value  of  DELMAX  is  only  used  in  cylindrical  finned 
analysis.  For  the  cylindrical  smooth  condenser,  DELMAX  is 
internally  generated. 
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APPENDIX  C 


SOURCE  CODE  LISTING 
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*****  SOLID  AND  FLUID  PROPERTIES  ****** 
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